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ABSTRACT 

Chiral perturbation theory in heavy-fermion formalism is developed for meson- 
exchange currents in nuclei and applied to nuclear axial-charge transitions. Cal- 
culation is performed to the next-to-leading order in chiral expansion which 
involves graphs up to one loop. The result turns out to be very simple. The 
previously conjectured notion of "chiral filter mechanism" in the time component 
of the nuclear axial current and the space component of the nuclear electromag- 
netic current is verified to that order. As a consequence, the phenomenologically 
observed soft-pion dominance in the nuclear process is given a simple interpre- 
tation in terms of chiral symmetry in nuclei. In this paper we focus on the 
axial current, relegating the electromagnetic current which can be treated in a 
similar way to a separate paper. We discuss the implication of our result on 
the enhanced axial-charge transitions observed in heavy nuclei and clarify the 
relationship between the phenomenological meson-exchange description and the 
chiral Lagrangian description. 
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1 Introduction 



By now there exist a large number of unambiguous experimental evidences M for 
meson-exchange currents in nuclear responses to electroweak probes. We also have avail- 
able a rather satisfactory and successful theory to describe the large bulk of experimental 
observations Q. While inherently phenomenological in character, the approaches taken so 
far to describe meson-exchange currents have been commensurate with the ingredients that 
account for our progressive understanding of nuclear forces and to the extent that nucleon- 
nucleon interactions are now fairly accurately understood, one can have a great deal of 
confidence in the theoretical tool with which the effect of exchange currents is calculated. 
There remains however the fundamental question as to how our phenomenological under- 
standing of nuclear forces and associated meson currents can be linked to the fundamental 
theory of strong interactions, QCD. 

In this paper we make a first step towards answering this question by applying chiral 
perturbation theory (ChPT) to nuclear electroweak processes. To start with, we assume 
that at low energies dominated by infrared properties of QCD, the most important aspect of 
QCD is the spontaneously broken chiral symmetry and hence that in nuclear dynamics, it is 
chiral symmetry that plays a predominant role Q . The important role of chiral symmetry 
in nuclear physics was recognized early on by Chemtob and Rho for exchange currents 
[Q] but this issue was recently given a stronger impetus and a more modern meaning by 
Weinberg in connection with nuclear forces || and by Rho || in connection with what is 
known "nuclear chiral filter phenomenon" (for definition, both intuitive and more rigorous, 
see later). 

The key question addressed here is this: To what extent can nuclear processes be 
described by QCD or equivalently at low energies by chiral perturbation theory? Weinberg 
approaches this issue by studying nuclear many-body forces. We propose here to do the 
same by looking at nuclear response functions responding to slowly varying electroweak 
fields 0. We suggest that chiral perturbation theory can be made - under certain conditions 
specified below - considerably more powerful and predictive for response functions than for 
nuclear forces. In calculating nuclear forces to loop orders in chiral perturbation theory, 
one encounters a plethora of counter terms to renormalize the theory, most of which are 
not accessible by experiments furthermore there are contact four-fermion interactions in 
the Lagrangian - most of which are again unknown parameters - that have to be carefully 
examined and treated. In principle, this may be feasible, perhaps with the help of lattice 
QCD calculations but in practice it may not be possible to make clear and useful predictions 
because of many uncontrolled parameters. While the tree order chiral theory justifies d 
posteriori the current nuclear physics practice of using two-body static forces |J, it appears 
that chiral symmetry will be unable to make any truly significant statement on the structure 
of nuclear forces for sometime to come. A major new development will be required before 
one can make a prediction that goes beyond the accuracy of the phenomenological approach 



1 



which has been strengthened by the wealth of experimental data. On the contrary, as we 
will show in this paper, when the formalism is applied to nuclear response functions, in 
particular, to exchange currents, it can make a highly nontrivial and potent prediction. 
This is because nuclear short-range correlations generated by nuclear interactions at short 
distance - which while poorly understood of their mechanism, are nonetheless operative 
in nuclear medium - screen all the contact interactions, both intrinsic and induced and 
consequently all of the four-fermion (and higher) counter term contributions, effectively 
"filtering" off the ill-understood short-range operators: Given phenomenological information 
on nuclear wave functions at short distance, the short-range suppression helps in simplifying 
nuclear response functions. In addition, in certain kinematic conditions, higher order chiral 
corrections are found to be naturally suppressed. The suppression of the many-fermion 
counter term contributions at the one-loop order that we are studying is, as will be stated 
more precisely later, a consequence of the fact that such terms occurring at high orders 
in the chiral expansion reflect the degrees of freedom that enter directly neither in nuclear 
forces nor in nuclear currents at the chiral order considered. The combination of these 



two phenomena lead to the "chiral filtering" proposed previously [10]. In this paper, we 
will establish this chiral filtering to one-loop order in chiral perturbation theory. This 
will provide, in our opinion, the very first compelling explanation of the pion-exchange 
dominance observed in axial charge transitions (considered in this paper) as well as in 
radiative np capture or in threshold electrodisintegration of the deuteron. 

As stressed by Weinberg, chiral perturbation theory is useful in nuclear physics only 
for "irreducible" diagrams that are by choice free of infrared divergences. This means that 
both in nuclear forces and in exchange currents, reducible graphs are to be taken care of by 
a Schrodinger equation or its relativistic generalization with the irreducible graphs entering 
as potentials. This also implies that in calculating exchange currents in ChPT, we are to use 
the wave functions so generated to calculate matrix elements to obtain physical amplitudes. 
This is of course the standard practice in the theory of meson-exchange currents but it is 
also in this sense that ChPT is predictive in nuclei. Clearly this precludes what one might 
call "fully consistent chiral perturbation theory" where nuclear forces, nuclear currents and 
wave functions are all calculated to the same order of chiral perturbation expansion. Such 
a calculation even if feasible is likely to make no sense. A little thought would persuade the 
reader that it is a futile exercise. 

We will here focus on the irreducible diagrams contributing to exchange currents. We 
will calculate next-to-leading order terms in the chiral counting involving one-loop graphs. 
In doing this, we will employ the recently developed heavy fermion formalism (HFF) [11]. 
The standard ChPT 11^] arranges terms in power of (d/A x ) and/or of (m n /A x ) where d is 
four-derivative acting on the Goldstone boson (viz, pion) field, m n the pion mass (~ 140 
MeV) and A x ~ 1 GeV, the chiral expansion scale. It has been established that this 
expansion works well at low energies for such processes as tttt scattering. However the 
situation is different when baryons are involved. The dynamically generated masses of the 
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baryons are of 0(A X ) and hence when the baryon field is acted upon by time-derivative, it 
gives an O(l) term. Therefore a straightforward derivative expansion fails. (Incidentally 
this is also the reason why a chiral Lagrangian describing pion interactions well with low- 
order derivative terms does not necessarily describe well skyrmion properties.) The HFF 
circumvents this difficulty in rearranging the derivative expansion. Indeed the principal 
virtue of the HFF is that it provides a consistent chiral expansion in Q/A x where Q is 
four-derivative on pion field or pion mass or space derivative on baryon field; it avoids time 
derivative on baryon field which is of order A x which is not small. The standard ChPT 
involving baryons jDJ can in principle be arranged to give a similar expansion. However 
it requires a laborious reshuffling of terms avoided in the HFF. The distinct advantage 
of the HFF is that the multitude of diagrams that appear in such calculations as ours in 
the standard ChPT involving baryons get reduced to a handful of manageable terms, thus 
alleviating markedly the labor involved. We will see that there is an enormous simplification 
in the number of terms and in their expressions. The potential disadvantage might be that 
the HFF is not fully justified for the mass corresponding to that of the nucleon and hence 
higher order "1/m" corrections may have to be systematically included. We will examine 
the class of approximations we make in the calculation by looking at the next order terms. 
It turns out that the leading "1/m" correction is absent in our calculation. We shall discuss 
this matter in the concluding section. 

While the procedure is practically the same, the resulting expression for electromag- 
netic (EM) current is somewhat more involved. We will therefore not treat it here although 
we shall give a general treatment of the theory applicable to both axial and EM currents. 
The detailed analysis on the EM currents, together with an application to threshold np 
capture, will be reported in a separate paper [14]. Both currents are intricately connected 
even at low energy through current algebras and we will need some vertices involving the 
EM current. 

It is perhaps obvious but we should stress that for both vector and axial-vector 
currents, relevant symmetries (i.e, conserved vector current and partially conserved axial- 
vector current) are preserved to the chiral order considered since both nuclear forces and 
currents are treated on the same footing with the same effective Lagrangian. More on this 
point later. 

This paper is organized as follows. In Section 2, we state our basic assumption in 
applying ChPT to nuclear dynamics. In Section 3, we describe the effective chiral La- 
grangian with which we develop heavy- fermion formalism including "1/m" corrections. We 
also define the relevant kinematics we will consider. The chiral counting rules are given in 
Section 4. In Section 5, the renormalization of n-point vertices that enter in the calculation 
is detailed. For the sake of making this paper as self-contained as possible and to define no- 
tations, we also list the renormalized quantities for the pion and the nucleon following from 
the Lagrangian. Readers familiar with renormalization of heavy-fermion chiral Lagrangian 
could proceed directly to subsection 5.4. Two-body exchange currents are calculated in 
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Section 6. Both momentum-space and coordinate-space formulas are given. Numerical 
analyses are described in Section 7. In Section 8, we explain why there are no other graphs 
that can contribute to the same order and point out in what circumstances they can show 
up in physical observables. Concluding remarks including those on the observed enhanced 
axial-charge transitions in heavy nuclei are made in Section 9. The Appendices A-I list all 
the formulas needed in the calculation. 

This paper is written in as a self-contained way as possible so as to be readable by 
those who are not familiar with the recent development in the field. Some of the material 
are quite standard and readily available in the literature. Most of them however serve as a 
check of our calculation. 



2 Strategies in Nuclear Physics 

We wish to calculate operators effective in nuclei for transitions induced by the vector 
and axial vector currents of electroweak interactions, denoted respectively by and 
associated with the electroweak fields and Au- In principle there will be n-body currents 
for ./V > n > 1 in iV-body systems. Here we will focus only on one- and two-body currents, 
ignoring those with n > 2. The reasons for so doing are given in the literature |2|, |J] but we 
will later show that n-body currents for n > 2 are suppressed to the order considered for 
long wavelength probes. 

The diagrams we wish to calculate are generically given by Fig. 1. They correspond 
to the standard definition of single-particle and two-particle exchange currents entering in 
the description of nuclear response functions to the external electroweak fields. These have 
been calculated before in terms of phenomenological Lagrangians. Here we wish to do so 
using chiral perturbation theory (ChPT), starting with a chiral Lagrangian that is supposed 
to model QCD at low energies. Following the chiral counting rule we will derive later, we will 
restrict our consideration to one-loop order, which corresponds to going to the next order 
in the chiral counting to the leading soft-pion limit. Although one-loop calculations have 
been done before for nucleon properties and for infinite nuclear matter [16], they have 



up to date not been performed in finite nuclear processes. We believe this work constitutes 
the first attempt to implement consistently chiral symmetry in nuclear processes. 

In dealing with divergences encountered in calculating loop graphs, in particular 
the loops involving two-pion exchange, we will need a certain prescription for handling 
operators that are short-ranged in coordinate space. This prescription does not follow from 
chiral symmetry alone and will have to be justified on a more general ground. Specifically, 
we argue that consistency with ChPT demands that zero-range interactions be "killed" by 
nuclear short-range correlations: 

• Firstly, the zero-range operators that come from finite counter terms appearing in 
four-fermion interactions figure neither directly nor indirectly - but importantly - in 
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the successful phenomenological nucleon-nucleon potentials and hence must represent 
the degrees of freedom unimportant for the length scale involved. In fact, one can 
show (see Appendix I) that the counter terms we need to introduce (denoted k\ ' 
later) cannot arise, unlike in the better understood it-it scattering jl?]], from an ap- 
proximation of taking an infinite mass limit of the strong interaction resonances such 
as the vector mesons p, u etc. which have a scale comparable to the chiral scale and 
play an important role in boson-exchange potential models. 

• Secondly, ChPT by its intrinsic limitation cannot possibly provide a nuclear force 
that can account for the interactions shorter-ranged than two pion or one vector- 
meson range at most. Thus the truly short-range interaction known to be present 
in the nucleon-nucleon interaction must involve elements that are not calculable by 
means of finite-order chiral expansion even if such an expansion existed. Thus it would 
be inconsistent to put a part of such interactions into the currents in the context of 
ChPT without a similar account in the nuclear force. It is known that even to one-loop 
order, the number of counter terms is so large in the calculation of nuclear forces that 
it is highly unlikely that one can make a meaningful prediction based strictly on low- 
order chiral perturbation expansion |8j . As suggested in Ref . || , one should implement 
ChPT calculations with phenomenological informations whenever available. 

• Applied to the "irreducible diagrams" that enter in the definition of exchange cur- 
rents, ChPT screens out the short-range part of the interaction which originates from 
dynamics of possibly non-chiral origin. When the matrix elements of the operators 
arising from the irreducible graphs are calculated with wave functions suitably com- 
puted in the presence of two-nucleon potentials, the short-range correlation built into 
nuclear wavefunctions must therefore suppress strongly interactions that occur at an 
internuclear distance < 0.6 fm, automatically "killing" the 5 function interactions 
associated with finite counter terms. This fact will be kept in mind when we derive 
two-body operators in coordinate space. 

We will present, at several places in the paper, arguments to justify the above procedure 
which purports to establish that the only unknown parameters in the theory must be (a) 
negligible in magnitude and (b) further suppressed by nuclear correlations when embedded 
in nuclear medium. 

There is nothing very much new in our calculation of the one-body operators except 
for its consistency with chiral invariance. As for the two-body exchange currents, our 
results are new. There are two graphs to consider: One-pion exchange (Fig. 2a) and two- 
pion exchange (Fig. 2b). Both involve one- loop order graphs. Note that we are to calculate 
only "irreducible graphs." 
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3 Effective Chiral Lagrangian 



We begin with the effective chiral Lagrangian that consists of pions and nucleons 
involving lowest derivative terms [JD]] relegating the role of other degrees of freedom such as 
vector mesons and nucleon resonances A to a later publication^ 1 , 



C = Nii-fidp + Tj-m + igA-fTsAjN -^C a (NT a N 
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+ ^Tr (v M StV^s) + ^M 2 F 2 Tr{E) + ■■■ + C CT , (1) 

where m ~ 939MeV is the nucleon mass, gA — 1-25 is the axial coupling constant and F ~ 93 
MeV is the pion decay constant. The ellipsis stands for higher derivative and/or symmetry- 
breaking terms which will given later as needed. We have written the Lagrangian with the 
renormalized parameters m, gA, F and M with suitable counter terms C CT to be specified 
later. Under the chiral SU(2)xSU(2) transformation^ 2 , the chiral field £ = exp(i^) 



transforms as £ —* g n T,g^ (<? R , <? L £ SU(2)) and the covariant derivative of the chiral field 
transforms as £ does, 

V M S = d^-i^ + A^ + i^-A^) 

-> gn^^gl (2) 

where the external gauge fields = • \ and A^ = A^ ■ \ transform locally 



V^ + A^^V^ + A'^ = <7r(V m + AJgt - id^g n ■ gt, 
V»-A^V'-A'= 5l (V m - A^gt - id^ ■ gt 



In our work, only the electroweak (SU(2) x U(l)) external fields will be considered. The 
Lagrangian of course has global SU (2) x SU (2) invariance in the absence of the pion mass 
term. Non-linear realization of chiral symmetry is expressed in terms of £ = \/E = exp(i^) 
and U = U(£,g h ,g n ) defined with £ 

e - 9n£rf = UCgt 



#1 While the vector mesons and the nucleon resonances (in particular, the A) play an important role in 
nuclear phenomenology - and they can be easily implemented in ChPT at least in low orders, they are 
unimportant for the process we discuss in this paper. It is not difficult to see which processes require such 
degrees of freedom but we will not pursue this matter, for a treatment of such processes goes beyond the 
framework of ChPT. 

#2 We are using a slightly unconventional notation of Ref. ju| which we will follow in this paper. This 
facilitates checking our results on single-nucleon properties against those derived in J]l| using standard 
(relativistic) chiral perturbation expansion. The more familiar transformation of the chiral field used in the 
literature is gotten by replacing E by E f . We are also working with the exponentiated (Sugawara) form of 
chiral Lagrangian instead of Weinberg's J5| used previously. They are of course equivalent. For the rest we 
will follow the Bjorken-Drell metric and convention. 
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Now nucleon field N transforms as N — > UN , and covariant derivatives of nucleon field and 



chiral field transform as nucleon field does, D^N — > UD^N and A^ — > C/A^C/f where ^ 



r 



" ~ 2 



A, = -^(V M S)^ = -p,94-V(V, + ^ + k(V,-i/ (3) 



Itt (y y\ tt _ I 

2^ v ^ / s 2 ^ 7 ~p^j 2^ v ' ^ 1 w ~p^ ■ 2' 

The U can be expressed as a complicated local function of £, g L and g R . The explicit form 
of U is not needed for our discussion. 

Note that we have included the four-fermion non-derivative contact term studied re- 
cently by Weinberg|5|. We will ignore possible four-fermion contact terms involving deriva- 
tives (except for counter terms encountered later) and quark mass terms since they are not 
relevant to the chiral order (in the sense defined precisely later) that we are working with. 
The explicit chiral symmetry breaking is included minimally in the form of the pion mass 
term. Higher order symmetry breaking terms do not play a role in our calculation. 

3.1 Heavy-fermion formalism 

For completeness - and to define our notations, we sketch here the basic element of the 
heavy-fermion formalism (HFF) fll8j| applied to nuclear systems as developed by Jenkins and 
Manohar[ll wherein the nucleon is treated as a heavy fermion. As stressed in Introduction, 



the relativistic formulation of ChPT works well when only mesons are involved but it does 
not work when baryons are involved since while space derivatives on baryons fields can be 
arranged to appear on the same footing as four-derivatives on pion fields, the time derivative 
on baryon fields picks up a term of order of the chiral symmetry breaking scale and hence 
cannot be used in the chiral counting. This problem is avoided in the HFF. To set up the 
HFF, the fermion momentum is written as 

= mv^ + W (4) 

where is the 4— velocity with v 2 = 1, and is the small residual momentum. (In the 
practical calculation that follows, we will choose the heavy-fermion rest frame = (1,0).) 
We define heavy fermion field B v (x) for a given four- velocity v^, by^ 5 



B v (x) =e imv - x N(x). (5) 
The field B v is divided into two parts which are eigenstates of tf, 

B v = fiW + = ^-^B v + ^--^B v = P+B v + P_ B v . (6) 



#3 We have defined two covariant derivatives involving chiral fields, V M E and A M . We can express one in 
terms of the other, but it is convenient as done frequently in the literature to use V M E for the meson sector 
and A M for the meson-nucleon sector. 

#4 Another familiar field redefinition is B v (x) = e iml "° v ' x N(x), with v 2 = 1. This definition gives exactly 
the same physics to the lowest order in — expansion. 
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As defined, B v +) (B ( v } ) can be identified as positive (negative) energy solution. As will be 
justified in the following subsection, the negative energy solution is suppressed for large 
baryon mass and its contribution can subsequently be incorporated as higher-order correc- 
tions in the inverse mass expansion. Thus to the leading order, the fermion loops can be 
ignored. With the neglect of the negative energy solutions, we have a useful relation for 
gamma matrices sandwiched between spinors which holds for any T, 



B V TB V = B v ffTB v = B V T^B V = B V ^{T, j6}B v . 
It follows from this identity that 

B vl5 B v = 0, B vl »B v = v^B v B v . 
Let us define spin operators S% by 

B vl5 ^B v = -2B V S%B V 

or explicitly 

The spin operators have the following identities, 



[S^S V V \ = ie^VaSp with e i23 = l- 



From the anti-commutation rule, we have 



S v ■ S v 



3 
~4 



(<? ■ s v ) 



SyS^S va — -(d — 3)S^ ~ -S 1 ^, 



(q-vf-q 2 



(7) 
(8) 

(9) 
(10) 

(11) 
(12) 

(13) 
(14) 
(15) 



where d is the dimension of the space-time, d = and we have defined e = (4 — d)/2. 
Between spinors, we have the approximate relations 



B, 



ByS^By ~ 

Dy, Oy 



Uy 



-<7 • v, -a 
2 ' 2 

i -, 
--v x a 
2 



(16) 
(17) 



with a the usual Pauli spin matrices. We are suppressed by a 

factor of v = O (Jji) where Q is the characteristic small momentum scale for processes with 
small three-velocity. Since 

[SS,£%] = l -(a^ + M^) 
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where a^ u = | [7^, j u ], we also have 

B V (^<T^B V = B V [SZ,S»]B V (18) 

B v (a^ 5 )B v = 2iB v (v»S»-v»S> J ')B v . (19) 

We are now in position to write down the chiral Lagrangian (||) in HFF. The nucleon 
part of the Lagrangian becomes 

N{i@-m)N = B v iv ■ dB v (20) 
and the corresponding nucleon propagator S(mv + k) is 



iS(mv + k) = — . (21) 

v ' vk + i0+ v ; 

Our chiral Lagrangian (|]) expressed in terms of the heavy-fermion field to leading (£e, 

zeroth) order in — takes the form 

1 /— \ 2 
-(" ' ~ " ~ ' 

2 



C = B v [iv fl (d fl + T^ + 2ig A S^A^B v --C a (B v T a B l 



F 2 1 
+ —TV (v^VE) + -M 2 F 2 TV(X!). (22) 

In practical calculations, the chiral field X or £ is expanded in power of the pion field. The 
explicit form resulting from such expansion as well as the vector and axial-vector currents 
calculated via Noether's theorem are given in Appendix A. 

3.2 — corrections 

m 

As mentioned, the HFF is based on simultaneous expansion in the chiral parameter 
and in "1/m". We have so far considered leading-order terms in 1/m, namely, 0((l/m) ). 
We now discuss — corrections following closely the discussion of Grinstein [|IJ| . We choose to 
do this in a perhaps more general way than needed for our purpose. Consider the following 
Lagrangian 

£ = N[ip-m + -fT B A ll ]N -^C a (NT a N) 2 



B [ijp _ m (i _ + 7 ^ 75 A M ] B - -C a {BT a Bf 



* 5 Although we do not actually encounter it in our calculation, it might be worthwhile to point out one 
technical subtlety. The infinitesimal i0 + is inserted to define the singularity structure of the propagator. 
When we encounter a d— dimensional loop integral we first perform the Wick-rotation to put it in the 
Euclidean metric. In doing this, we assume that the first and third quadrants (in the plane of real 1° vs. 
imaginary 1°) contain no poles. If we can take the flow direction of the loop-momentum to be the direction 
of the fermion momentum, there is no problem. However for some graphs it is impossible to do this. For 
instance consider a two-nucleon box diagram. In this case, we have one fermion line in which the loop- 
momentum flow direction is opposite to that of the fermion arrow. In this case, the fermion propagator is 
of the form 

= r + 2iTvS(v ■ k). 

v ■ k - i0+ v ■ k + i0+ y ' 
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where we have included Weinberg's four-fermion contact term with T a an arbitrary her- 
mitian operator which we assume to contain no derivatives. We have also introduced an 
arbitrary "axial" field which we take to be hermitian and free of gamma matrices. (Here 
and in what follows, we shall omit the subscript v in B v , B v and S%.) The equation of 
motion satisfied by B is 

[G - m(l -f)]B = (23) 



with 

G = g-C a r a (BT a B), g = + -f^. 
Multiplying (^) on the left by P- , we obtain 

P-GB- 2mB { - ) = 



(24) 



which leads to 



Now multiplying P + to (|23[), we get 



= — P- GB = — P_GB {+) + O (W . 
2m 2m 



m- 



(25) 



P+GB = 



which gives 



G+ — GP-G 
2m 



(26) 



Given this, it is now a simple matter to write down the Lagrangian that gives rise to the 
equation of motion to the desired order The result correct to 0(l/m) is 



C 



B 



9 + 77-5 p - 9 
2m 



B 



1 



r a + 77- (raP-g + gP-r a ] 

2m 



B 



with g defined in To put this into a more standard form, we use the identities 



(27) 



P+gP+ 

P + gP-g'P+ 
P + gP-gP + 



P+-y,g} P+, 

P + {P-,g}g'P + = P + g{P^,g'} P 4 



Pa 



\{^g 2 }-\{g,gy,g}} 



Pa-. 



(28) 



One can show from these identities that 

(b {+) r a bW) ® p (+) r a p_ =o (29) 

* 6 We remind the reader that one should not insert the solution of into the original Lagrangian, since 
in Lagrangian approach, what is important is the form, not the value. However in Hamiltonian approach, 
the insertion of the solution for into the original Hamiltonian is allowed. 
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for any T a = {1, 75, 7^, 7^75, 07^}. This allows us to simplify the Lagrangian further to 
the form 

C = B{iv-D + 2S-A)B-\c a (BT a B) 2 + ^-B(-D 2 + {v-D) 2 



2 " V " J 2m 
+ [S», S u ] [Dp, D u ] - (v ■ A) 2 -2%{vA,S- D}) B, (30) 

for general T a allowed by symmetries. In our case, A^ = igA^^ = it/AgC (V^E)^, so our 
^ term Lagrangian is 

SC = -^B(-D 1 + {ii-Df + \S^,S"]\D IL ,D t ,\+ a \{v &f + 2g A {v A, S • £>}) B 

+ o(-jr)- < 31 ) 



While Eq.(31) is the first "1/m" correction, it is not the entire story to the order 
considered. One can see that it is also the next order in the chiral counting in derivatives 
and it is expected in any case, independently of the inverse baryon mass corrections. Thus 
in a practical sense, the coefficients that appear in each term are parameters rather than 
fixed by chiral symmetry in HFF. We also note that in the derivation given above, neither 
0{l/m) corrections to the quartic fermion term nor sixth-order fermion terms are generated 
by the (1/m) expansion. This of course does not mean that such terms cannot contribute 
on a general ground. We have however checked that no such terms arise from exchanges of 
single heavy mesons that are formally "integrated out" which means that our calculations 
are not affected by such terms, at least to the order we are concerned with. 



4 Counting Rules 

In this section, we rederive and generalize somewhat Weinberg's counting rule|| 
using HFF. Although we do not consider explicitly the vector-meson degrees of freedom, 
we include them here in addition to pions and nucleons. Much of what we obtain turn out 
to be valid in the presence of vector mesons. Now in dealing with them, their masses - 
which are comparable to the chiral scale A x - will be regarded as heavy compared to the 
momentum probed Q - say, scale of external three momenta or m^. 

In establishing the counting rule, we make the following key assumptions: Every 
intermediate meson (whether heavy or light) carries a four-momentum of order of Q. In 
addition we assume that for any loop, the effective cut-off in the loop integration is of order 
of Q. We will be more precise as to what this means physically when we discuss specific 
processes, for this clarifies the limitation of the chiral expansion scheme. 

An arbitrary Feynman graph can be characterized by the number En (Eh) of external 
- both incoming and outgoing - nucleon (vector-meson) lines, the number L of loops, the 
number 1^(1^, Ih) of internal nucleon (pion, vector- meson) lines. Each vertex can in turn 
be characterized by the number di of derivatives and/or of m n factors and the number rij 
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(hi) of nucleon (vector-meson) lines attached to the vertex. Now for a nucleon intermediate 
state of momentum = m/ + where k^ = 0(Q), we acquire a factor Q _1 since 

S F (mv + k) = ^— = 0{Q- 1 ). (32) 
v ■ k 

An internal pion line contributes a factor Q~ 2 since 

A^;ml) = -^ = 0(Q^) (33) 

while a vector- meson intermediate state contributes Q° (~ O(l)) as one can see from its 
propagator 

Af(? 2 ; m 2 v ) = 1 2 ~ — L- = O(Q ) (34) 
<^ — m y — my 

where my represents a generic mass of vector mesons. Finally a loop contributes a factor Q 4 
because its effective cut-off is assumed to be of order of Q. We thus arrive at the counting 
rule that an arbitrary graph is characterized by the factor Q" with 

u = -I N - 2I n + 4L + Y,di (35) 

i 

where the sum over i runs over all vertices of the graph. Using the identities, I n + Ih + In = 
L + V — 1, Iff = \ J2i hi — 4r and Tjv = \ J2i n % ~ ^2~> we can rewrite the counting rule 

E N + 2E H ^ _ m + 2hi 
u = 2 \-2L + 2_^Vi, Vi = di-\ 2. (36) 

i 

We recover the counting rule derived by Weinberg || if we set Eh = hi = 0. 

The situation is different depending upon whether or not there is external gauge field 
(i.e., electroweak field) present in the process. In its absence (as in nuclear forces), Vi is 
non-negative 

Hi + 2hj , . 

di + ■ 2 - 2 > 0. (37) 

This is guaranteed by chiral symmetry |5|. This means that the leading order effect comes 
from graphs with vertices satisfying 

m + 2hj , , 

di + 1 2 * -2 = 0. (38) 

Examples of vertices of this kind are: 7r k NN with k > 1 (di = 1, rii = 2, hi = 0), 
ZiiViV (di = 0, raj = 2, hi = 1), (ATW) 2 (d t = 0, m = 4, ^ = 0), hir k with fc > 1 (dj = 
1, Hi = 0, /ij = 1), etc where h denotes vector-meson fields. 

In iV./V scattering or in nuclear forces, En = 4 and Eh = 0, and so we have v > 0. The 
leading order contribution corresponds to v = 0, coming from three classes of diagrams; 
one-pion-exchange, one-vector-meson-exchange and four-fermion contact graphs. In ttN 
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scattering, En = 2 and Eh = 0, we have v > 1 and the leading order comes from nucleon 
Born graphs, seagull graphs and one-vector-meson-exchange graphs.^ 7 
In the presence of external fields, the condition becomes || 

4 + !*±^-2) >-l. (39) 

The difference from the previous case comes from the fact that a derivative is replaced 
by a gauge field. The equality holds only when hi = 0, in = 2 or hi = 0, n; L = 0. We 
will later show that this is related to the "chiral filter" phenomenon. The condition ( |39"| ) 
plays an important role in determining exchange currents. Apart from the usual nucleon 
Born terms which are in the class of "reducible" graphs and hence do not enter into our 
consideration, we have two graphs that contribute in the leading order to the exchange 
current: the "seagull" graphs and "pion-pole" graphs both of which involve a vertex 



with Vi = —1. On the other hand, a vector-meson-exchange graph involves a ^ = +1 
vertex. This is because di = 1, hi = 2 at the J^hh vertex. Therefore vector-exchange 
graphs are suppressed by power of Q 2 . This counting rule is the basis for establishing the 
chiral filtering even when vector mesons are present (see Appendix I). Thus the results we 
obtain without explicit vector mesons are valid more generally. 

5 Renormalization in Heavy-Fermion Formalism 

In this section, we discuss renormalization in heavy fermion formalism. Most of the 
renormalized quantities that we will write down here have been obtained by others in stan- 
dard ChPT [13]. We rederive them for completeness and as a check of our renormalization 



procedure. 

For reasons stated above, fermion loops are suppressed in HFF. Our basic premise 
is that antiparticle solutions should be irrelevant to physical processes in large-mass and 
low-energy situations. Their effects can however be systematically taken into account in 
"1/m" expansion. 

o o o o 

We shall denote "bare" quantities by m, M, F, 9a and the corresponding "renor- 
malized (to their physical values)" quantities by m = m + 5m, M, F, qa, respectively, 
for nucleon mass, pion mass, pion decay constant (~ 93 MeV) and axial coupling constant 
(~ 1.25). 

5.1 Dimensional regularization 

We adopt the dimensional regularization scheme to handle ultraviolet singularities in 
our loop calculations. It has the advantage of avoiding power divergences like 5(0) ~ A^ ut 

* 7 We note here that scalar glueball fields \ play only a minor role in tyN scattering because the y7T7r vertex 
(di = 2, rii = 0, hi — 1) acquires an additional Q power. 
#8 These are standard jargons in the literature. See 
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where A cu t is the cut-off mass. In d = 4 — 2e dimensions, all the infinities are absorbed in -. 
When heavy fields are involved, somewhat different parametrization schemes and integral 
formulas are needed. The relevant ones for our calculation are 



A n B r(n) Jo {2XA + BY+ 1 



A n BC r(n) 7o 7o + zB + (1 - z)C] n+2 



and 



oo 



r(n) / dAA fc (A 2 + a 2 y n = -a^-^T ( n - ^) V ( . (42) 



o 



This integral is singular when 



„-*±l = 0,-1,-2,..- (43) 



so e must be kept finite until the integration in A is performed. Some relevant integral 
identities needed in this paper are given in Appendix C. 

As is customary in the dimensional regularization, we introduce an arbitrary mass 
scale fj,. After renormalization, the results should, of course, be independent of the scale \i. 
Here some comments are in order regarding the one-loop renormalization scheme. First, all 
the divergences of our theory can be classified in two classes by their degrees of divergence: 
quadratic and logarithmic divergences. The quadratic divergences are removed by counter 
terms that are of exactly the same form as the lowest chiral- order Lagrangian, Co as given by 
(|2~2|). Thus these quadratic singularities can be absorbed into the renormalization process 
of the basic quantities, namely, gA, F, M, m. To take care of the logarithmic singularities, 
we include counter terms that are higher chiral order by Q 2 than Cq. As it is an arduous 
task |l3j to write down all possible counter terms, we shall write down only the counter 
terms needed for the calculation. All the quadratic divergences can be written in terms of 
A(M 2 ) and all the logarithmic divergences in terms of r] defined by 

a / n jri\ A d f d d l 1 M 2 .[M 2 \~ e 

A(M 2 ) = — -.- =- r(-l + e w 1 



i J (2ir) d l 2 -M 2 16tt 2 v ' ^vr/i 2 

M 2 ( 1 



16vr 2 



€ 



+ l + r'(l)-ln— )+0(e), (44) 




1 . . / M 2 \ 1/1 . M 2 \ . 

i = i^ r(£ »(4^j = i^l7 +r(1) - ll, wj + 0(£) (45) 

where e = (4 — d)/2 and r'(l) ~ —0.577215. Note that both A and r\ are scale-dependent 
and singular. But the coefficients of the counter terms (written down below) are also scale- 
dependent and singular. To remove the scale dependence and singularity, one must adjust 
the coefficients of the counter terms. This procedure however is not unique since finite 
parts of the coefficients of the counter terms are totally arbitrary. There are many ways to 
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eliminate this non-uniqueness leading to a variety of "subtraction schemes." In this paper, 
we use the scheme whereby renormalization is made at the on-shell point for the nucleon 
and at zero four- momentum for the pion and the current. Thus the quantities g A , F and 
the coefficients of the counter terms are defined at zero momentum of the axial current and 
the pion. 

To make our discussion on renormalization streamlined, we list below all the counter 
terms needed in our work (for the meson sector, see @) to which we will refer back as 
required; 



+ 



B 



-5m + (Z N - l)iv -D + ij^(v Df 



B 



Co 



ig A ^BS^v a [ A^D u D a - D v \ L D a + D u D a A M 



B 



I — ttBv 



F 2 ' 



D^,[D a ,Dp] 



B 



{vD,[D a ,Dp}}B 



+ l jk 

+ ij^B[v-A,[vD,vA]]B 
+ i^B[S-A,[vD,S-A]]B 

9A 4 1} Br a D„B ■ B [v ■ A, r ] S»B 



8F 4 

_ SA J2) 
4i?4 

+ h.c.} 



df (B [5 m , S v ) D^B ■ Bv ■ A S U B + B [S^, S u ]v A D^B ■ BS V B) 



(46) 



with 



5m 

Cl 

C3 
C 4 
Q> 
C6 



1 + 



3(d - l)g\ A(M S 
4 F 2 



32vrF 2 ' 

\aAJl + cf , 
d- 3 



3 

1 r 
6 



9 2 av + <£, 



l + (d+l)<& 



v + 4, 



-2g 2 AV + cf, 
V + 4, 



-{d-S)g%7i + ($ 
(d-l)g 



.(2) 



8g 2 A v 



(47) 



where cf (i = 1,2, ■■-,6) and k^~' 2 ^ are finite renormalized constants that we will refer 
to as "finite counter terms." Since these finite counter terms and the finite parts of loop 
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contributions are scale- independent, our final results also are scale-independent and regular. 
The chiral counting is immediate from the counter-term Lagrangian. We should mention 

(1 2) 

for later purpose that the two-derivative four-fermion counter terms proportional to d\ 
cannot be gotten from single low-mass resonance exchanges and hence do not figure in 
long-range as well as intermediate-range NN potentials. 

We should note here that although the above counter-term Lagrangian contains the 
isospin matrix r, chiral invariance of the Lagrangian is preserved. This can be verified by 
noting that 

(R'^aR).. {R~ l r«R) kl = fa) fi (r«) w (48) 

where (i,j,k,l = 1, 2) and R is the SU(2) matrix for chiral transformation, B — > fl-B.p 9 ! Also 
note that is hermitian while A^ and [S' /i ,5 1 '] are antihermitian and 7°-D / ! t 7 =D^=d^ 
— where is the antihermitian operator defined by = <9 M + T^. 



5.2 Pion properties to one loop 



Since due to pair suppression fermion loops can be ignored, renormalization in the 
pion properties is rather simple. The wavefunction renormalization Z n , renormalized mass 
M and pion decay constant F (here as well as in what follows renormalized at q 2 = with 
M^O) are given by PlP°l 

2 A(M 2 ) 

3 F 2 ' 



1 



M 



o 2 

M 



F 



1 + 



_ A(M 2 
" ~~2F 2 

A(M 2 )' 
F 2 



(49) 
(50) 

(51) 



with A(M 2 ) defined in Eq.(@). 



5.3 Nucleon properties to one loop 

One-loop graphs for nucleon propagator are given in Fig. 3. Fig. 36 vanishes due to 
isospin symmetry, so only Fig. 3a survives to contribute to the nucleon self-energy E, 

E(^ = -3||^/ Mi + ^_ M2) (52) 



where (and in what follows) 



d d l . . 

(53) 



(2?r) d ' 

#9 Actually this equation can be simply understood if one uses an 0(3) representation. In this case, R is 
an orthogonal real matrix and (r a )ij = —itaij, a,i,j = 1,2,3. 

* 10 We have not put in the counter terms that appear at 0(Q 4 ) in the pure meson sector, so the L\ terms of 
Gasser and Leutwyler jjjj are missing from our expression. Since the meson sector proper does not concern 
us here, we will leave out such terms from now on. 
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From this, we get 

3(d- l)g\ A(M 2 ) , , 

Z N = 1 + E'(0) = 1 + V 4 ;yA y p2 ' (54) 

where the prime on E stands for derivative with respect to u ■ fc. In our case confined to 
irreducible graphs, there is no nucleon pole, so we can set v ■ k = in the denominator of 
the nucleon propagator. For an off-shell nucleon, we have 

E(vk) = (d-l)^h(vk) (55) 

with the function h(v ■ k) defined by 
f 11 

I v(l + k)(l 2 -M 2 ) = 9apKv • k) + VaM ' • (56) 

where (• • •) stands for a function that does not concern us. The evaluation of h(v ■ k) is 
described in Appendix D. For (v ■ k) 2 < M 2 , we have 

m = 3{d ^ mm 2 ) y - y 3 + M (M 2 _ y2) %)i (57) 



where y = v ■ k and rj is a singular quantity given by ( |45| ) and /io is a finite function, the 
explicit form of which is given in Appendix D. We note that the above equation contains 
the - divergence in the coefficient of (v ■ k) 3 as well as in the coefficient of (v ■ k) (i.e, in 
A(M 2 )). This additional singularity arises also in the conventional method. See the paper 



by Gasser et. a/.[13|. The counter term needed to remove this divergence, as given in (|46|), 



is 



S CT (y) = 5m-(Z N -l)y + ^y 3 . (58) 
In order to regularize the propagator subject to the condition E(0) = E'(0) = 0, we choose 

3(d-l)g 2 A A(M 2 ) 



Z N = 1 + 
5m 



4 F 2 

3g\M 3 



The result is 



32ttF 2 ' 
3 



••5/r 



E(y) = 6m + ± y 3 + Jf (m 2 - y 2 ) h (y). (59) 

Here the finite constant cf is in principle to be determined from experiments. To see its 
physical meaning, we should look at a process involving an off-shell nucleon. For instance, 
when v ■ k = ±M, we have 

Ft 

E(±M) = ±%M 3 + 5m. (60) 
F z 

Finally, the ^ correction is readily seen to be 

SS(k) = — — \k 2 - (v ■ k) 2 ] . (61) 
2m L J 
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5.4 Renormalization of 3-Point Vertex Functions 



In this subsection, we shall calculate three-point vertex functions to one-loop order, 
in particular, J^NN and irNN given in Fig. 4, where J M = A^iy^) denotes the axial vector 
(vector) current. We treat the vector current simultaneously since some vertices involving 
it figure in our calculation. Each vertex function is a sum of contributions from tree graphs, 
one-loop graphs, wavefunction renormalization, higher-order counter term insertion and 
O (j^j corrections, if needed. The tree-graph contribution to J^NN is 

r ANN( tree ) = 9AT a S^, 

r^(tree) = ^ (62) 

Wavefunction renormalization produces multiplicative coefficients Z, i.e., Zn for T A ^ N anc ^ 
zh N for T a nNN etc. 

Unless noted otherwise we will always set the momentum flow of all pions and currents 
to be outgoing. The current-off-shell-nucleon couplings that we will consider are of the type 

N(mv + k) -> N(mv + k - q) + J%(q), 
N(mv + k) — > N(mv + k — q) + ir a (q) 

with the relevant momenta indicated in the parentheses. 

For completeness, we list in Appendix F all the contributions to the three-point 
functions of each Feynman graph. 

5.4.1 Axial vertex function T A ^ N 
For off-shell nucleons, we find 



~^ANN ~ 9 A T a S^ 1 



A(M 2 ) , 3(d-l)gi 2 d-3g\ 

1 + p2 + ^ A ( M ) + —4~p2 h 3( V ■ k > V ■ «) 



(63) 



where hs(v ■ k, v ■ q) = ^- [h(v ■ k — v ■ q) — h(v ■ k)] is evaluated in Appendix D. The singu- 
larity in the above equation is removed by the counter term contribution 

[r^Ar] CT = 9 A r a S» {-^) [3(« • k f - 3v ■ kv ■ q + (v ■ q) 2 ] (64) 

with 

c 2 = ~^-9aV + 4 (65) 

where c R is a finite renormalized coupling constant. Adding the counter term contribution 
to the loop contribution, we obtain a renormalized axial coupling constant by g^ = 9a{^ = 
0, q = 0) where gA(k, q) is defined by 

r^N R (k,Q) = 9A(k,q)T a S^. (66) 
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Physically gA(k,q) is just the axial charge form factor for the incoming nucleon of momen- 
tum mv^ + k^ and the axial current carrying the momentum q^. Explicitly it is given 

by 

^^ = l + {^Mv-k,v-q)-^[3(v-k) 2 -3v-kv-q + (v-q) 2 } (67) 
where h 3 (v ■ k, v • q) is a finite function defined by 



h 3 (v ■ k,v -q) = h 3 (v ■ k,v ■ q) + A(M 2 ) - -rj 3(v • k) 2 - 3v • k v • q + (v • q) 2 



and 



9 A = 9a 



1 + 



A(M 2 

'f 2 ~ 



~ 1.25 



(68) 



(69) 



where we have equated the renormalized g A to the experimental value. Finally, the ^ 
correction is 

ST^nn = -^9AT a v»(2k-q)-S. (70) 

For on-shell nucleons, v ■ k = v ■ q = or in Breit frame, v ■ k = \v ■ q, so we find 

g A {k,q)=g A . (71) 

Note that we have neither momentum dependence nor — corrections for an on-shell nucleon 
or in Breit frame. This means that there will be no one-loop correction to the irNN vertex 
in the exchange currents calculated below. 

Given an experimental axial charge form factor of the off-shell nucleon, one can fix 
the constant cf. The vertex T% NN can be obtained by a direct calculation or by means of 
Lehmann-Symanzik-Zimmermann (LSZ) formulation: Both give the same result 

■ 9A{k,q) 



-pa R 
1 ivNN 



(k,q) 



T a q- S. 



(72) 



5.4.2 Vector vertex function Ty^ N 

The full form of Ty^^ is rather involved, 



where 



1 VNN 



1 

" T 2 
(d-i) 

+2 



f F ' + h [q ' s > s " ] F2 ~ q " F *) 



h(q 2 )-A(M 2 



+ 



A 

F 2 



^A(M 2 )-h 3 (vk,vq)-B 2 (k,q) 



(v • q) 2 - q 2 



d 1 

^B 2 (k,q)+AvqB 1 (k,q)j , 



FY = 



Am^B (k,q), 
v ' q f („ 2 \ + & ' 

-prMQ ) + ^2 



(d+l)-2((v.q) 2 -q 2 ) 



_d_ 

8q 2 



B 1 (k,q) 



(73) 



(74) 
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with the functions fi, hi and F>i given explicitly in Appendix B, D and E, respectively. 
Although the above equations appear to have quadratic divergences, they actually have 
only logarithmic divergences as can be seen below : 

,2 



F 



1 



_ (2 „. *_„.,)!!_+_!_, 



2 d + 1 gi 2 
' + —z-rSiMT 



6 F 2 



+ - 



d- 



1 9 A 



3(v ■ k) 2 — 3v ■ k v ■ q + (v ■ q)" 



+ 



M q\ q 2 A , 



\ + {d+l)g\ 



v ■ q 
6F 2 



rj + • 



(75) 



where the ellipsis denotes finite and 0(Q n )\ n> 2 terms. Quadratic divergences disappear 



0) 



because of EM gauge invariance. We see that Ty'^ N (k = 
logarithmic divergences are removed by the counter term in Eq. ( |46| ) of the form 



l 

2 'a 



T a v^. The remaining 



I 1 VNNIct 



_C3_ 

2F 2 
2F 2 



3(vkf 



■•o 



,2^ 



3v ■ q v ■ k + (v ■ 



T a v(2k-q) [q-S,S»] 



with 



(76) 



C4 



l + {d+l)g A r, + <% 



-2g\ V + cf. 



R 



(77) 



When v ■ k = v ■ q = 0, we have -F3 



and 

J2 



16vr 2 F 2 



2(l + 2<&)ir 2 (g 2 



*2V) 



g\ m 

'AttF 2 



j dz^M 2 - z(l - z)q 2 + 1 



(78) 
(79) 



where we have added the correction appearing in the second term of F^ , It is easy 
to see in (|78|) that the counter term constant c§ can be related directly to the isovector 
charge radius of the nucleon. We will give the precise relation later. Ki(q 2 ) are the finite 
pieces of the functions fi(q 2 ) defined in Appendix B, 

q 2 ~ 

dz In 1 — z(l — z) 



K (q 2 ) 
K 2 (q 2 ) 



dz z(l — z) In 



1 



M 2 
z(l 



z) W 2 



(80) 



* 11 Although the loop contribution to the Pauli form factor is finite and hence requires no infinite 
counter term, there is a finite counter term contributing to it which we did not - but should - include in our 
formula. As pointed out by Bernard et al. |2(J, the finite counter term can be considered as coming from 
the p exchange as in vector dominance picture. This point will be addressed more precisely in |mJ. 
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5.5 Renormalization of 4-Point Vertex Function 

In this section, we study the 4-point vertex functions denoted and T^y b to one 
loop as given in Fig. 5, corresponding to the process 

N(mv + k)^ N(mv + k-q a -q b ) + J£{q a ) + n b (q b ) 

where the isospin components of the current and the pion are denoted by the subscripts a 
and b, respectively. Here v ■ k represents how much off-shell the incoming nucleon is and 
v ■ (k — q a — q b ) the same for the outgoing nucleon. For tree graphs, we have 



r^f(tree) = -^e abc T c v^. 



(81) 



The full formulas for non-vanishing graphs for off-shell nucleons are given in Appendix 
G. Here we limit ourselves only to the on-shell nucleon case. For axial-charge transitions, 
only the six graphs Fig. 5 (a) — (/) survive. Figures 5 (g),(h) are proportional to Su, so 
suppressed for the time component and Figures 5 (i) — (n) are proportional to v ■ S = 0. 
Figures 5 (o) — (r) do not contribute to the axial current. (We have included them for later 
purpose, see [0].) 

Adding the loop contributions and tree graphs with their wavefunction renormaliza- 
tion constant, we have 

^nA = e abcT c ^ A + i^ab 

(82) 



with 



1 nA 



~ {l - ^ [(l + (d - l)g\) h{q 2 ) - 8g 2 A (q ■ S) 2 f 2 (q 2 ) - v q a h£(v ■ q a ) 
2 A [q . S ,S»]B {q 2 ) 
+ H (^T- Ub ■ S, S»} hf(v ■q a ) + 6 ^ [q b ■ S, S»] hf(v ■ q a )) , (83) 

+ 2(q.S) 2 



p3 y 4 imo j '"4 \- ^aj i 4 

r 1 — d _ „ „ s o d > (- 

B (q ) + -j^v ■ q a h {v 



K'l = -(g a + 3g 6 y^ 



3F 3 



4 ^ ' dq 



-3^ {q b • S, S»} hf(v ■q a ) + d -^ [q b ■ S, S»] hf(v ■ (84) 

where q = q a + q b , f\{q 2 ) = fi{q 2 ) - /i(0) and h^ ,A (y) = \ [hi(y) ± hi(-y)]. The integrals 
defining the functions fi(q 2 ) for i = 0, 1,2,3 are listed and evaluated in Appendix B. (/ij's 
are defined in Appendix D and -Bo m Appendix E.) The log divergences contained in these 
vertices are removed by the counter term contribution, 

1 



1 nA 



-t-abc p \FvNn1ct 
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i c 5 «/ \2 



with 



C5 = 77 + c f> 



c 6 = -(d-3)g\ri + cg (86) 

where are renormalized finite constants listed in Eq.(|46|) and [ry-'^ rAr ] CT is given by (|76|) . 

With soft momentum, we have v ■ q a = for which we obtain a surprisingly simple 
expression, viz, 

K=A' ab = -tabcT^F? \t) (87) 

where t = q 2 = (q a + q^) 2 and F± is the isovector Dirac form factor Eq.(|78]). The one-loop 
renormalization of the irANN vertex corresponds to the isovector charge radius obtainable 
from the form factor F± for which the finite counter term c§ plays a key role. We see that 
and T^nn are re l a ted, respectively, to Ty^ N and T^y b calculated in Appendix G. That 
the A^ttNN vertex for a soft pion is simply given by FY has of course been understood 
since a long time via current algebra and also in terms of the /9-meson exchange. 

Finally for the — correction, one can readily obtain the corrections to the vertices 

*r£f = -e abc r c ^- F (2k^-q>'-v^2vk-vq) + 2[q-S,S^) 
-iS ab v^v -q a ' J 1 



AmF 

where k is the residual momentum of the incoming nucleon and k — q = k — q a — q b is that of 
the outgoing nucleon. An important point to note from Eq.(p8|) is that for the case v-q a = 
and on-shell nucleons, we have no contribution from — corrections to the time component 
of the axial current and the space component of the vector current. This is the basis for the 
pair suppression we will exploit in the application to axial-charge transitions in nuclei. 

The complete list of the four-point functions involving the vector and axial-vector 



currents needed here and in [14] is given in Appendix G. 



6 Two-Body Exchange Currents 

So far we have computed one-loop corrections to the graphs involving one nucleon. 
They are extractable from experimental data on nucleon properties. In this section we 
incorporate the above corrections into - and derive - two-body exchange currents in heavy- 
fermion chiral perturbation theory. As shown previously J/J, the time component of the 
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axial- vector current (and also the space component of the vector current [0]) in the long- 
wavelength limit is best amenable to chiral perturbation loop calculations. ^ 1 ^ We will work 
out the computation to one loop order corresponding to the next-to-leading order in the 
chiral counting rules as derived in Section 4. The process of interest is 

N( Pl ) + N( P2 ) - N(p[) + N(p' 2 ) + J£(fc), 

where we have indicated the relevant kinematics with q 2 = p 2 — P2, Qi = Vi ~ Pi an d the 
energy-momentum conservation q\ + q 2 + k = 0. The process is soft in the sense that 

(Q 2 \ 

v ■ k ~ v ■ qi = v ■ (p'i — pi) ~ O — ~ 

\m I 

where m is the nucleon mass or chiral scale. This kinematics markedly simplifies the calcu- 
lation. Clearly this kinematics does not hold, say, for energetic real photons. 

It is convenient to classify graphs by the current vertex involved. The graphs that 
contain J^irNN play a dominant role since J^irNN ~ for the axial current (and 
J^irNN ~ S 1 ^ for the vector current). The graphs which contain J^NN (and J^ttitNN 
for the vector current) can be ignored to the relevant order because the role of the vector 
and axial currents is interchanged. In what follows, we discuss the axial current only. The 



argument for the vector current goes almost in parallel and will be detailed in [14]. 

What we are particularly interested in is the time component of the axial current, 
with axial-charge transitions in nuclei in mind. This is where the "pion dominance" is 
particularly cleanly exhibited. The space component is also interesting both theoretically 
and phenomenologically. Theoretically Gamow- Teller transitions - observed to be quenched 
- represent the other side of the coin relating to the chiral filter phenomenon discussed above, 
namely that chiral symmetry alone or more precisely soft mechanisms associated with it 



cannot make statements on this quantity [21]. Empirically the quenching phenomenon is 



closely associated with the missing strength of giant Gamow- Teller resonances in nuclei. 
Since the treatment of the space part of the axial current requires going beyond chiral 
perturbation theory, we will not pursue this issue any further in this paper. 
For convenience, we define an "axial-charge" operator A4 by^ 13 



#12 This is the crucial point in using ChPT in nuclei that quantifies the general discussion given in Section 
2. Since this point is often misunderstood by nuclear physics community, we would like to stress it once 
more although to some it may sound obvious and repetitive. The chiral counting on which our analysis is 
based is meaningful only if « 1 where Q is the characteristic momentum or energy scale involved in the 
process. Therefore ChPT cannot describe processes that involve energy or momentum scale exceeding that 
criterion. This means that processes probing short-distance interactions are not accessible by finite-order 
ChPT. In particular two-body currents describing an internuclear distance ri2 < 0.6 fm cannot be probed 
by the expansion we are using. We argued before - and will make use of the fact - that there is a natural 
cut-off provided by short-range nuclear interactions that go beyond the strategy of ChPT which allows a 
meaningful use of the small Q expansion. 

* 13 The operator Ad is an isovector but in what follows we will not explicitly carry the isospin index. 
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We decompose M. into Ai = Mtree + Mioop where Mtree denotes the axial charge operator 
coming from the one-pion-exchange tree graph (i.e, the soft-pion term) and Mi oop is what 
comes from loop corrections. Further we decompose Mi oop into 

Mioop = M ln + M 27 r (90) 

where Mi n denotes the loop correction to the one-pion-exchange axial charge operator 
(also referred in the literature to as "seagull graph") and M 2 -k the contribution from two- 
pion-exchange graphs and tree graphs involving four-fermion contact terms with counter- 
term insertions. We will later argue that the latter does not contribute. One-loop graphs 
involving four-fermion contact interactions, while allowed in the relevant chiral order, do 
not however contribute either. 



6.1 Results in momentum space 

As stated above, the non-zero contributions to the time component of the axial cur- 
rent come only from the graphs that contain a J^ttNN vertex. The tree seagull graph 
supplemented with a vertex form factor - and properly renormalized - leads to^^ 



A? ree + Al = in x r 2 ^ Q2 • S 2 -^—^FYiql) + (1 - 2) (91) 

where F± is the Dirac isovector form factor of the nucleon. To the order considered, there are 
no further corrections. The present formalism allows us to calculate within the scheme the 
form factor FY ■ The tree-graph (or soft-pion) contribution corresponds to taking FY = 1. 
The difference (FY — 1) is given by Figs. 5(a) — (/). Taking — ► 0, we encounter two 
spin-isospin operators, 

= -2in XT 2 q-Sx + (1^2)~iT 1 X r 2 q- {S x + <? 2 ), (92) 
= 2(fi+f 2 ) [g-S2,Si-S2] + (1^2)~i(fi + f 2 )<f-<TiX«T 2 (93) 

with qV- = q% — — q% . With the help of these operators, we can rewrite ( |9l| ) as 

Mk = _T(1) M^? [*?V)-l] (94) 

With Mtree = tJ—,. 



As for two-pion-exchange and four-fermion contact interaction contributions, the rel- 
evant diagrams are those given in Fig. 6(a) — (k) and their symmetrized ones. Before going 

#14 We denote particle indices by i — 1, 2 without expliciting heavy fermion fields. For instance, Si should be 
understood as the spin operator sandwiched between B v and B v of particle 1 with velocity v. In this section, 
q 2 is the four-momentum squared of the pion but we are concerned with the situation where << 1, so 
the static approximation q 2 ~ — |q| 2 will be made in practical calculations and also in Fourier-transforming 
to coordinate space later. In fact, the static approximation is not only natural for the chiral counting but 
also essential for suppression of n-body forces and currents for n > 2. More on this later. 
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into any details, one can readily see that each graph in Fig. 6 contributes a term of at least 
0(Q). This can be shown both in the conventional method and in HFF by observing that 
their contributions vanish if we set M = qff = k^ = 0. This assures us that our counting 
rule is indeed correct. Therefore we can neglect all the graphs proportional to since the 
axial-charge operator involves S° ~ 0(Q/mjv) as stated before. Figures (/), (g), (h) and 
(j) belong to this class. Now Fig. (e) is identically zero because of the isospin symmetry 
and Fig. (i) is proportional to v ■ S = 0. The graph (k), involving time-ordered pion 
propagators, are the so-called "recoil graphs" j| which as we shall argue in Section 8 will 
be cancelled by similar recoil terms in reducible graphs. So we are left with only the four 
graphs (a), (6), (c) and (d) to calculate. Without any further approximation than using 
HFF, the full expression of the four graphs comes out to be 



A" (a) 
A»(b) 

A»(d) 



-(2r 2 

(2r 2 + iri x f 2 
(-2ri 



-* \ QA 

IT1 X T 2 ) ^ 
9A 



8F 4 
2r 2 + if\ x r 2 



{v»qz-S\ f (q z 2 ), 
(q2-Sv»\ f (q 2 2 ), 



3 

(2T! + 2f 2 + ir x x f 2 )^ (S a v»\ (S"S' 3 ) 2 I^(q 2 



with fo(q 2 ) given in detail in Appendix B and I^^piq) defined by 



I n,ap(q) 



hvlv (l + q)(l 2 - M 2 ) [{I + q) 2 - M 2 } ' 
This integral is evaluated in Appendix E. Using the conditions v ■ qi = and k^ 
can rewrite them in a symmetrized form 



A^ia + b) 
A^(c + d) 

A" 



9A 



16vr 2 8F 4 



K (q 2 



16ir 2 i] 



rd), 



(95) 

(96) 
0, we 

(97) 



{[-(d- 1)16tt 2 77 + 3K (q 2 ) + 2K 1 (q 2 )] T« 



16tt 2 16F 4 
-8 [K (q 2 ) - 167^7?] T( 2 )}, 
.m 9A ( d (D T (D + d (2) T ( 2 ) 



where q^ = i(g 2 — q\Y and 



16F 4 



At" 



(98) 



is the contribution from the counter-term Lagrangian 



The Ki(q 2 ) are finite pieces of the functions fi{q 2 ) defined in Appendix B, i.e. 



K (q 2 



K,{q 2 



dz In 



dz- 



z{1 ~ z) w 



-z(l 



z)q 



' M 2 - z(l - z)q 2 ' 

The expressions (]97| ) and ((9^) contain singularities in r], which are removed by the counter 
term contribution with 



d 



(i) 

4 
(2) 



.(2) 



+ 



(d-l)9. 



(99) 
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where the renormalized constants K4's are finite and scale- independent. 

The resulting two-body axial-charge operator including finite counter-term contribu- 
tions is 



M 2 * = 



2) 



16n 2 F 2 

4 ^(^ (1) + 4 2) t (2) )- (ioo) 



The two-body axial-charge operator due to loop correction is then the sum of (|9J) and ( fOO ) 



Mioop = M l7T + M 27r . (101) 

As it stands, the constants K4's are the only unknowns in the theory as they cannot be 
determined from nucleon-nucleon interactions as mentioned before. They could in principle 
be extracted from two-nucleon processes like N + N — > N+N+ir but they appear as higher- 
order corrections and it is inconceivable to obtain an information on these presumably small 

(i) 

constants from such processes. However as argued above, we expect the constants k\ to 
be numerically small and what is more significant, when we go to coordinate space as we 
shall do below to apply the operator to finite nuclei, they become 5 functions and will be 
completely suppressed as we discussed in Section 2. In momentum space, such constant 
terms have also to be removed as done for the celebrated Lorentz-Lorenz effect (or more 
generally for the Landau-Migdal g' ) in pion- nuclear scattering |^| . It should be stressed 
that once the constant counter terms are removed, no unknown parameters enter at next to 
the leading order in the chiral expansion in nuclei. It is also noteworthy that to the order 
considered, the loop contributions are renormalization-scale independent. 



6.2 Going to coordinate space 

Applications in nuclear transitions are made more readily in configuration space. 
Furthermore considerations based on ranges of nucleon-nucleon interactions which seem 
necessary for rendering chiral symmetry meaningful in nuclei are more transparent in this 
space. Therefore we wish to Fourier-transform the operators (|97| ) and ( |98[) into a form 
suitable for calculations with realistic nuclear wave functions. In doing this, we will treat 
the pion propagator in static approximation, namely, q 2 ~ — |q| 2 . In Appendix B, we show 
how the highly oscillating integrals involved in the calculation can be converted into inte- 
grals of smooth functions by performing the Fourier transform before doing the parametric 
integration. Since the spin-isospin operators T^ and T^) contain q — which is a derivative 
operator in configuration space, it is convenient to define 

T (1) = f\ x f 2 f • (a i + cr 2 ), 

= (fi + r 2 ) f ■ {ax x a 2 ). (102) 
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Writing Eqs.(|94|) and ( |100|) in coordinate space which we will denote by Ai to distinguish 
from the momentum-space expression, we obtain - modulo 5 function terms mentioned 
above - the principal result of this paper: 



Mtree{r) 
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dr 
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4nr 
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(103) 

iT (r) - K (r)] + (2 + 4<&) [# 2 ( r ) _ £ 2 ( r )] j , 

(104) 

T«+2< 7 iK (r)T( 2 ) 1,(105) 



3fl^ - 2 . . 1 „ 



VWi T (r) +7W 2 ^(r). 



(106) 



As defined, A4 n7T are n7r exchange corrections to the soft-pion (tree) term. Mi 00 p is therefore 
the total loop correction we wish to calculate. The explicit forms of the functions Ki(r) 
and Ki(r) are given in Appendix B. 

As noted above, the constant c§ can be extracted from the isovector Dirac form factor 
of the nucleon, i.e., 



R M 2 M 2 , 9 , v 



0.04784. (107) 
It is interesting to separate what we might call "long wavelength contribution" from M\ n , 
Min(r) = 5 so ftMtree(r) + (short range part) 

where 



Osoft 



M 2 



° 3 F 2 + 16vr 2 F 2 



1 + 3g 2 



7T 



_ _(l + 2ff l) -- — 



17 



7T 



9 



~ 0.051 (108) 



and compare this one-loop prediction for S so ft * 15 to what one would expect from the 
phenomenological dipole form factor 

2 



fY(^ 



A 2 



A 2 



(109) 



with A = 840 MeV. This form factor leads to the following one-pion exchange contribution 
to M-ioap-, corresponding to (104): 

2 

I -1 



p~i dipole 



fa) 



M 2 



A 2 



A 2 -M 2 



* 15 It is worth noting that this contribution is generic in the sense that it is more or less model-independent: 
It is of the same form and magnitude whether it is given by chiral one-loop graphs or by the vector dominance 
(see Appendix I) or by the phenomenological dipole form factor. This may have to do with the fact that 
it is controlled entirely by chiral symmetry. It is curious though that this longest wavelength effect is an 
enhancement rather than quenching usually associated with form factors. 
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1 



A 2 



2 \ A 2 - M 2 



A 2 + 



A 2 



A 2 -M 2 



1 A\ 
r z r J 



• (no) 



Identifying the first term of ( |11C| ) with the first term of ( |104j ), we see that 8 so ft corresponds 
to 

( A 2 \ 2 



A 2 - M 2 



1 ~ 0.0571. 



It is remarkable - and pleasing - that the one-loop calculation of S so ft is so close to the 
empirical value. Furthermore the remaining term in ( |104| ) involving the functions K[ cor- 
responds - and when applied to the process of interest, is numerically close - to the second 
(short-ranged) term in (110). 



7 Numerical Results 



In order to get a qualitative idea of the size involved, let us first look at the magnitude 
of the relevant terms given in momentum space. For this purpose we set q 2 ~ — |q| 2 ~ — Q 2 , 
where Q is taken to be a characteristic small momentum scale probed in the process which we 
take to be of order of at most. For convenience we shall factor out the tree contribution 



from the expression ( 101 ) and write it as 

M =M t ree(l+6 M + 0(Q 3 )) 



fill) 



where 5 M is the chiral correction of 0(Q 2 ) that we have computed (relative to the tree 
contribution). We obtain 



where, setting (T^) = (T^) in nuclear medium and dropping the K4's 

Q 2 



J2w 



4F 2 

_Q 2 + M 2 
4F 2- 



2 fV) V + I±^ (Q 2 ) 



l + 2gl 
2vr 2 



K 2 (Q 2 



5gj + 2 
16vr 2 



K (Q 2 



For Q ~ m7r ~ 140 MeV, g A = 1.25 and F ss 93 MeV, we get 

\&Xk\ ~ 0.045, |5 2 7r| ~ 7.5 x 10" 3 . 



(112) 



(113) 



This is consistent with the notion that at the relevant scale Q, the chiral correction remains 
small. 



* 16 One can show in fermi-gas model, Wigner's SU (4) supermultiplet model or even jj-coupling shell model 
of nucleus with one particle outside of closed core, (T^ 1 ') ~ (T^). This relation will be assumed in all 
numerical calculations that follow. We would like to thank Kuniharu Kubodera for his help on this relation. 

#17 Dropping the constant terms in momentum space is not fully justified unless all other terms of the same 
nature are removed as well. This problem is avoided in coordinate space. We give only the absolute values 
for the 8 m for the same reason. See below for more on this matter. 
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We now turn to a more realistic estimate of the chiral correction appropriate to the 
actual situation in finite nuclei. Calculating nuclear transition matrix elements in momen- 
tum space is cumbersome and delicate. There are several reasons for this. The most serious 
problem is the implementation of the short-range correlation. In the well-studied case as 
in the 7r-nuclear scattering, we know how to proceed, obtaining the celebrated Lorentz- 
Lorenz effect. Roughly the argument goes as follows | |22| |. Consider a term of the form 
q 2 /(q 2 + M 2 ) that figures in the p-wave pion-nuclear scattering amplitude, or more specif- 
ically in the interaction between the particle-hole states excited by the pion. Rewrite this 
as 1 — M 2 /{q 2 + M 2 ). Removing the constant 1 corresponds to suppressing a 5 function 
in coordinate space and leads to the Lorentz-Lorenz factor. Note that this procedure of 
accounting for short-range correlations can even change the sign. Unfortunately our case 
does not lend itself to a simple treatment of this kind because of the nonanalytic terms 
coming from the loop contributions: there is no economical way of "removing 5 functions" 
from them. This task is much simpler and more straightforward in coordinate space. 



Let us therefore turn to the coordinate space operators (|103| ) and ( 106 ). In Fig. 7, we 
plot Aitree (P-03|) and Aiioop ( |106|) as function of the internuclear distance r = \ri— rz\ setting 
= T^ 2 ) = 1. Some of the important features discussed in the preceding sections can 
be seen in this plot. While negligible at large distance, say, r > 1 fm, the loop corrections 
get progressively significant at shorter distances and at r ~ 0.4 fm, they are comparable to 
the soft-pion result. There is nothing surprising or disturbing about this feature at short 
distances. At shorter distances which are probed by the momentum scale approaching 
the chiral scale, there is no reason to ignore the degrees of freedom integrated out from 
the theory. Low-order calculations with higher chiral-order degrees of freedom eliminated 
cannot possibly describe the short-distance physics properly. This may be construed as a 
sign that ChPT is not predictive in nuclei. We claim that this is not so. The point is that 
as long as the scale Q probed by experiments is much less than the chiral scale, truncating 
higher chiral-order and shorter wavelength degrees of freedom as done in ChPT can be 
meaningful provided short-range nuclear correlations are implemented in the way discussed 
above. 

Calculations of the nuclear matrix elements with sophisticated wave functions in 
finite (light and heavy) nuclei - and comparison with experimental data - will be made 
and reported in a separate paper. Here for our purpose of getting a semi-quantitative 
idea, the fermi-gas model as used by Delorme |2^] will suffice. One could incorporate 
accurate correlation functions - and this will be done for specific transitions in finite nuclei. 
Here we will not do so. We shall instead take the simplest correlation function, namely 



g(r, d) = 9(r — d) with the cut-off distance d ~ 0.7 fm as used by Towner [24]. Since this is 
a rather crude approach, we will consider the range of d values between 0.5 and 0.7 fm. 

Specifically we are interested in the ratio of the matrix elements (M.i oop ) / '{M-tree} 
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which in fermi-gas model takes the form (see Appendix H) 

R (d n) = {Mloop) = ( 1U ) 

(Mtree) / d °° dr V [jl(»r)] 2 M tr ee(r) 

where pp and p = -^zF§ are, respectively, the fermi-momentum and density of the sys- 
tem, ji(x) = - and Mi oop {r) = Mi n (r) + M2n(r). Note that w(pp,r) = 
4irr [ji(per)] 2 /j§ can be viewed as a weighting function. Since this calculation is straightfor- 
ward, we shall not go into details here. For completeness, however, we sketch the calculation 
in Appendix H. 

In Fig. 8 are plotted the functions w{p ei r)M. tree {r) and w(p F ,r)Mi oop (r) with T"W = 
T( 2 ) = 1 for pp ~ 1.36 fm™ 1 corresponding to nuclear matter density. The ratio R(d, p) is 
plotted in Figure 9 for d = 0.5, 0.7 fm. For d = 0.7 fm which was used by Towner J24J, 
the loop correction is at most of the order of 10% of the soft-pion term at nuclear matter 
density. There are two important points to note in the result. The first is that separately 
the loop corrections to the one-pion term (i.e, M\ n ) and the two-pion term (i.e, M.2-k) can 
be substantial but the sum is small. The second point is that the resulting loop contribution 
has a remarkably weak density dependence. The first is a consequence of chiral symmetry 
reminiscent of the tree-order cancellation in linear a model of the nucleon pair term and 
the a-exchange term in the S-wave ttN scattering amplitude. The second observation has a 
significant ramification on the mass dependence of axial-charge transitions in heavy nuclei 
to which we will return shortly. 



8 Other Contributions 

Here we briefly discuss what other graphs could potentially contribute and the reason 
why they are suppressed in our calculation. Consider the two-body graphs given by Figure 
6(A;) where the pion propagators are time-ordered. They belong to what one calls "recoil 
graphs" in the literature 0. To 0(Q 2 ) relative to the soft-pion term, these graphs - and 
more generally all recoil graphs including one-pion exchange - do not contribute. The reason 
is identical to the suppression of three-body forces as discussed by WeinbergQ: the graphs 
in Fig. 10 are exactly cancelled by the recoil corrections to the iterated one-pion exchange 
graphs that are included in the class of reducible graphs. Thus to the extent that the static 
approximation is used in defining the one-pion exchange potential, these graphs should not 
be included as corrections. Incidentally this justifies the standard practice of ignoring recoil 
graphs in calculating exchange contributions in both weak and electromagnetic processes 
in nuclei. 

We have ignored in our calculation three-body and higher-body contributions such 
as Figure 11. The reason for ignoring these graphs is identical to that used for proving the 
suppression of three-body and other multi-body forces ||. As in nuclear forces, they can 
contribute at 0(Q 3 ) relative to the soft-pion termQ. 
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An interesting question to ask is in what situations the approximations that justify 
dropping the graphs considered here break down in nuclei. It is clear that the static ap- 
proximation - one of the essential ingredients of the heavy-fermion formalism - must break 
down when the energy transfer involved is large. Imagine that one is exciting a A resonance 
in nuclei by electroweak field. The energy transfer is of the order of 300 MeV, so the static 
approximation for the pion propagator involving a 7rAiV vertex cannot be valid. In such 
cases, one would expect that multi-body forces and currents suppressed in this work could 
become important. This suggests specifically that in electron scattering from nuclei with 
sufficiently large energy transfer, n-body currents (for n > 2) will become progressively 
more important in heavier nuclei. Combined with the dropping mass effect (i.e, "Brown- 
Rho" scaling mentioned below), one expects a large deviation from the standard mean- field 
description used currently. 

9 Conclusions and Discussions 

We have used heavy-baryon chiral perturbation formalism to calculate the leading 
corrections to the soft-pion axial-charge operator in nuclei. Exploiting short-range suppres- 
sion of the counter terms and other short-range components of the two-body operator, we 
have shown that the chiral filter mechanism holds in nuclear matter with a possible uncer- 
tainty of no more than 10%, thus confirming the dominance of the soft-pion exchange. In 
a separate paper, we will show that the same holds in electromagnetic responses in nuclei. 
Since the currents (both vector and axial-vector) are calculated consistently with the symme- 
tries involved, they are fully consistent with nuclear forces that are calculated to the same 
chiral order. Ward-Takahashi identities will be formally satisfied although in practice ap- 
proximations made for calculations may disturb them. The final consistency will of course 
have to be checked a posteriori case-by-case. 

Taking this result as a statement of chiral symmetry of QCD in nuclei, what can one 
learn from this concerning the phenomenological models popular in nuclear physics where 
one uses exchanges of all the low-lying bosons in fitting nucleon-nucleon scattering (such 
as the Bonn potential) as well as calculating the exchange currents? Suppose we denote 
the axial-charge two-body operator from one-pion exchange with form factors by A^, one- 
heavy- meson exchange with form factors by Ah, the axial current form factor by App, all 
calculated within a phenomenological model, then our result implies that for the model to 
be consistent with chiral symmetry, then the total must sum to 

A tatal =A ln + A H + A FF + ---^A soft (l + d), \S\«1 (115) 

where A so f t is the soft-pion term as defined in this paper and S is the next-to-leading term 
of 0(Q 2 ). Our calculation illustrates how individually significant terms conspire to give a 
small 0(Q 2 ) correction which is insensitive to nuclear density. 
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One other outcome of our result is that while a subset of graphs can have a substantial 
density dependence, the small net chiral correction from the totality of the graphs does 
not have an appreciable nuclear density dependence, at least in fermi-gas model. We see 
no reason why this weak density dependence should not persist in more realistic nuclear 
models. Thus assuming that n-th order chiral corrections for n > 3 (relative to the leading 
soft-pion term) are not anomalously large, we come to the conclusion that meson-exchange 
axial-charge contributions to nuclear matrix elements cannot be substantially enhanced in 
heavy nuclei over that in light nuclei. The question arises then as to what could be the 
explanation for Warburton's recent observation that while the mesonic effect is about 50% 
in light nuclei, it is required to be 100% in heavy nuclei such as in lead region ||25[|. One 



suggestion [26| was that the parameters of the basic chiral Lagrangian have to be modified 
in the presence of nuclear matter consistent with trace anomaly of QCD [^7|. It predicted 
that hadron masses and pion decay constants that appear in the single-particle and one- 
pion exchange two-body operators are scaled by a universal factor $ that depends on matter 
density. Another suggestion [24, 28] was that exchanges of heavy mesons a, p, uj, a\ etc could 
become important in heavy nuclei while relatively unimportant in light systems. The latter 
mechanism relied on nucleon-antinucleon pair terms in phenomenological Lagrangians. Both 
mechanisms seemed to qualitatively account for the enhancement. 

We wish to understand the possible link, if any, between the chiral Lagrangian ap- 
proach and the phenomenological approach that includes pair terms involving heavy mesons. 
Since within the chiral approach developed in this paper the pair is naturally suppressed 
as required by chiral symmetry and multi-body currents are also suppressed as discussed 
above, the scaling mass effect of [26, ^7]] is the only plausible mechanism left within low- 
order chiral expansion for the medium enhancement noted by Warburton. Needless to say, 
we cannot rule out - though we deem highly unlikely - the possibility that higher order chi- 
ral terms supply the needed density effect. Incorporating the possible 10% loop correction 
calculated above in the two-body operator and the scaling factor $ = m* N /miy 0.8, one 
gets in the scheme of Ref . [ p6j [ the enhancement in heavy nuclei (at nuclear matter density) 
e MEC ~ 2.1 which is reasonable in the lead region compared with the experimental value 



2.01 ± 0.05 [25]. Within the scheme, this is the entire story and a surprisingly simple one. 
Of course more detailed finite nuclei calculations will be needed to make a truly meaningful 
test of the theory. 



In the phenomenological approach studied by the authors in [24] and [28], there is no 
fundamental reason to suppress NN pairs, so that heavy mesons could contribute through 
the pair term. However the exchange of heavy mesons, particularly that of vector mesons, is 
suppressed by short-range correlations in nuclear wave functions. Furthermore in the model 



of Towner [24], a large cancellation takes place in the sum such that the relation ( |115| ) seems 
to hold weil p9[] : Towner finds 5 < 10% over a wide range of nuclei. It is naturally tempting 
to suggest that Towner's model gives a result close to ours because it is consistent with chiral 
symmetry, at least to the same order of chiral expansion as ours, with higher-order terms 
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implicit in Towner's model which need not be consistent with ChPT somehow cancelling 
out^ 15 . There is however one aspect that needs to be clarified: in the models of [24] and [28], 



there is a pair term associated with a scalar meson (a) exchange. In the chiral Lagrangian 
used in this paper, there is no equivalent scalar field. We have however the scalar field x 



associated with the trace anomaly of QCD which plays a role in the Brown-Rho scaling [27]. 
We believe that these two effects are roughly related in the sense discussed in Ref. [30]. 
In this sense we would say that the pair term involving the a meson is simulating the 
density-dependence of the nucleon mass in the one-body axial-charge operator. There is no 



mechanism in Ref. [24, 28], however, for the density dependence of Ref. [p6| in the two-body 
operator. We suggest that this can be generated by taking three-body terms with an NN 
pair coupled to a cr-exchange. 

An obvious omission in our treatment of the axial current is the space component of 
the current governing Gamow- Teller transitions in nuclei (and the time part of the electro- 
magnetic current in Jl4|). The reason for this was already stated at several points in the 
paper: this part of the current is not dominated by a soft-pion exchange and indeed as noted 



many years ago 21 ] it is rather the very short-ranged part of nuclear interactions (roughly 
equivalent to the removal of the 5 functions associated with the counter terms in the spin- 
isospin channel) that plays an important role, e.g. in quenching the axial-vector coupling 
constant from the free-space value qa = 1-25 to g* A ~ 1 in nuclear matter. (For a simi- 
lar situation with the isoscalar axial-charge transition mediated by a neutral weak current 
where soft-pion exchange is forbidden, see Ref.|31|].) Furthermore, three-body operators for 
the space component of the axial current may not be negligible. For instance, as one can 
see from Appendix A, the three-body Gamow- Teller operator involving one nucleon with an 
A^ttttNN vertex with the pions absorbed by two other nucleons is not trivially suppressed 
as it is for the axial-charge operator. This suggests that low-order chiral perturbation the- 
ory may have little to say about this aspect of nuclear interactions. It is intriguing that in 
nuclei, both chiral and non-chiral aspects of QCD seem to coexist in the same low-energy 
domain. This makes QCD in nuclei quite different from and considerably more intricate 
than QCD in elementary particles studied by particle physicists. 

Finally we mention a few additional issues we have not treated in this paper but we 
consider to be important topics for future studies. 

• It would be interesting to see what two-loop (and hopefully higher-loop) and corre- 
sponding chiral corrections do to the chiral filter phenomenon. Two-loop calculations 
are in general a horrendous task but the situation in nuclear axial-charge transitions 



* 18 The following observation may be relevant to our argument that the counter terms must be ignored. 
Suppose one constructs a purely phenomenological theory based on meson-exchange picture by fitting ex- 
periments but conform to the symmetries of the strong interactions. Towner's model is one such example. 
One can convince oneself that in such a model, it is not possible to generate counter terms of the K4 type in 
infinite mass limit. Therefore if such terms existed, then they must be due to degrees of freedom that are 
not relevant at the accuracy required. 
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might be considerably simpler than in other processes. 



It would be important to see whether ChPT is predictive for processes involving larger 
momentum transfers as well as large energy transfer. From our experience with the 
electrodisintegration of the deuteron at large momentum transfers where the naive 
soft-pion approximation seems still to work fairly well, we conjecture that the chiral 
filter mechanism holds still in some channels even in processes probed at shorter 
distances or at larger momentum transfers. But as mentioned above, large energy 
transfer electron scattering might require multi-body currents in heavy nuclei. 

In this paper, we worked with an effective Lagrangian in which all other degrees of 
freedom than pions and nucleons have been integrated out. It would be important to 
reformulate ChPT using a Lagrangian that contains vector mesons incorporated a la 
hidden gauge symmetry (HGS) |}2] fitJ-9 and also nucleon resonances (such as A). As 



mentioned before, we believe that the chiral filter argument presented in this paper 
is not modified in the presence of these resonances in the Lagrangian. In Appendix 
I, we show that the presence of vector mesons does not modify our prediction on the 
axial-charge operator. Furthermore it is not difficult to see that the baryon resonances 
- in particular the A resonances - do not contribute to the axial-charge transitions 
to the order considered. However as is known for Gamow- Teller transitions in nuclei 



[21], certain processes in nuclei might require, even at zero momentum transfer, an 



explicit role of some of these heavier particles. As recently shown by Harada and 



Yamawaki [34], vector mesons introduced via HGS can easily be quantized, so their 



implementation in ChPT would pose no great difficulty. 

A systematic higher-loop chiral perturbation approach using the same heavy-fermion 
formalism to kaon-nuclear interactions and kaon condensation has not yet been worked 
out. This is an important issue for hypernuclear physics, relativistic heavy-ion physics 
and stellar collapse [35]. 



Finally if the parameters of effective chiral Lagrangians scale as a function of nuclear 



matter density as suggested by Brown and Rho[p7|, then one expects that as matter 
density increases, many-body currents will become increasingly important even at 
small energy transfer. This was already noticed in [26] where the soft-pion exchange 



charge operator became stronger in heavier nuclei. We already noted that this effect 
will show up more prominently in nuclear electromagnetic responses with large energy 
transfer. Future accurate experiments in electron scattering off nuclei will test this 
prediction. Of course this issue has to be treated together with many-body forces that 
enter into such processes. 



* 19 As stressed by Georgi, the HGS is an approach most suited to a systematic chiral counting when vector 
mesons are explicitly present. See [p3[]. 
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Appendix A: Chiral Lagrangian Eq.([22|) Expanded 

For completeness we expand Cq in powers of the pion field and in external fields. We will 
group Cq by the number of external gauge field, Cq = Cq + Cq + Cq, 



Ml 



C 1 



C 2 



+ Biv ■ dB 



I 



6F 2 



vf 2 a ( vf • d^vf - (vr • d„7r) 2 
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M 2 
4!i?2 



(*■)= 



+ B 
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From the (A.2), we extract Noether currents, 
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Appendix B: Functions /i(^ 2 ) 



The functions fi(q 2 ) (i = 0, 1, 2, 3) fi guring in subsections (5.4.2) and (5.5) are defined 

by 

HQ 2 ) = ' ' 



{I 2 - M 2 ) [(I + g) 2 - M 2 ] ! 



/ 3 (9 2 ) = 2/ 2 (g 2 ) + i/ ( (? 2 ) (B.f) 
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where we have defined f^{q 2 ) through 

la (2/ + q)fi 



(I 2 - M 2 ) [{I + q) 2 - M 2 } 



g a pfi{q 2 ) + q a qph{q 2 )- 



(B.2) 



Here and in what follows, the mass M could be thought of as the pion mass M. One can 
verify that 

dfi(q 2 " 



dq 2 



f2(q Z )- 



After performing the parametric integration, we have 

1 



fo(q 2 



V 



16vr 2 



h{q 2 ) = A(M 2 )-V 



f2(q 2 ) 



6" 16vr 2 



2A'2(<f 2 ) - ^ofa 2 ) 



71 1 
4 + 



6 16vr 2 



where 
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and K(q )'s are finite functions given explicitly by 
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(B.4) 
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We should note that all the functions given above are positive definite for negative q 2 and 
vanish when q 2 = 0. For —q 2 <^ M 2 , we have 



K (q 2 ) = K x {q 2 



6 M 2 



+ -he), 



M 4 



M 4 
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In the chiral limit (—q 2 ^> M 2 ), they simplify to 



K (q 2 ) 
Ki(q 2 ) 
K 2 (q 2 ) 



In 

1, 
1 



M 2 



2. 



5 

18' 



In order to go to the r— space, we must Fourier-transform Ki(q 2 ) and M 2 j{M 2 — 
q 2 )Ki(q 2 ). This involves an integration^ 20 of highly oscillating functions. Instead of in- 
troducing a regulating function which kills contributions from large g 2 and performing a 
tricky numerical calculation, we transform the problem into an integral of a smooth function 
with the use of residue calculation. The point is that we do the parametric integration (of 
variable z) at the last step. To see how this work, let us rewrite Ki(q 2 ), 



K (-Q 2 ) 
K^-Q 2 ) 

K 2 (-Q 2 ) 



dx In I 1 + 



E 2 



dx 



Q 2 



Q 2 + E 2 ' 



dx — 



(B.8) 



where Q = |q| = \/—q 2 and E = E{x) = 2M/\/l — x 2 .^ 1 First we Fourier-transform 



algebraically the integrands of the above integrals and then do the parametric integral 
numerically. The Fourier transform of K\ becomes an elementary residue calculation with 
a pole at Q = iM, 

KAt) = S(t) - -!- / dxE 2 e~ Er (B.9) 
47rr Jo 

while Kq and K 2 are somewhat involved due to the logarithmic function. We rewrite the 
logarithmic functions (of Kq and K 2 ) into a simple pole form by integration by part 
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With the above equation and 
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we obtain the expressions for K(r, A) defined by 
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#20 Recall that Ki(q 2 ) does not go to zero when q goes to infinity. 

* ;21 We have made the change of variable, x — 2z — 1, to render the expressions more symmetric. 
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K n (r A) = — Cdx^- — (A 2 e~ Ar - E 2 e~ Er 

01 ' j 4^r J d 1 - x 2 A 2 - E 2 V A e ^ 6 

tf 2 (r,A) = -i- fdx-^ 9 fAV Ar -£V & ). (B.13) 

v 7 4vrr i 1 - x 2 ^4 12 J A 2 - £ 2 V J v 7 

Here the parameter A is introduced to regularize the integrals near the origin of the config- 
uration space. When A = oo, we get the expressions for the Ki(r) needed in the paper, 

KAr)= lim KAr, A). (B.14) 

A— >oo 

The above integrals are non-singular even near x = 1, since E increases so as to make the 
integrals regular. However the expressions contain highly singular terms near the origin of 
configuration space when A goes to infinity. To see this, note 

Ko (r) = --!- f 1 dx E 2 e~ Er + lim -L A 2 K (q 2 = A 2 ) e" Af \ 



4irr Jo 1 — x 2 A^oo Airr 

and 

lim — A 2 e" Ar = 5(r). 
A^oo Airr 

Now ii"o(A 2 ) goes to infinity logarithmically when A 2 goes to infinity. So, roughly, the 
second term of ( p2|) behaves in the limit of infinite A as 

A 2 

lim In— — 7<5(r). 

The mathematical reason for this behavior is not hard to see : For large Q 2 , both Ko and 
K2 increase logarithmically with a gentle slope. Thus they can be viewed as a constant at 
large Q 2 with their value increasing to infinity. The constant behavior leads to the S(r) 
function. This term singular at the origin is of course "killed" by the short-range correlation 
present in nuclear wave functions. 

In terms of the function so defined, we can immediately obtain Ki(r) defined by 

(4 e,q ' rifi( " Q W = K ^ M) - (B - 15) 

In summary we have the expressions for Ki(r) (valid for r > 0) 

**> - -SFjf* -is)*-"* <B16) 

and the expressions for Ki(r) (valid in the whole region) 
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Appendix C: Integral Identities 



In this Appendix, we list some useful identities for the integrals we need to evaluate. 
Consider the following integral without imposing the condition v 2 = 1, 



I a (v,M 2 ) = f - 
Ji '< 



lr 



vi{i 2 -M 2 y 

then Lorentz covariance implies that its most general form is 

I a (v,M 2 ) = v a I (v 2 ,M 2 ). 
Now multiplying v a to both sides, we get 

A(M 2 ) =v 2 I (v 2 ,M 2 ) 



or 
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livl{l 2 -M 2 ) v 2 ' 

Using this expression, we obtain the following identities by successive differentiation with 
respect to v, 
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and so on. 



Appendix D: Functions h(v • /c)'s 

In this Appendix, we define - and give explicit forms of - the function h(v ■ k) that 
figures in the self-energy for the off-shell nucleon discussed in section (4.3). The basic one 
is ho(v • k), 
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(D.2) 



where y = v • k, = \J M 2 — y 2 , M = \J y 2 — M 2 and — ^ < sin _1 (x) < ^. We observe 
that the imaginary part appears only when v ■ k > M. The even part of ho has a very 
simple form 

ho(y) + h (-y) 



mv) = 



(D.3) 
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For special values of y, we have h Q (0) = -ff, h' (0) = 2rj, h (M) = -h (-M) = j^rj- 
The finite function ho(y) is defined by 

h{y) = Hv) ~ 2 vv- (D.4) 

We now examine the function h(v ■ k) defined by 
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If we multiply the above equation by g a ^ and vr , we obtain the following identity 



Kv) 



1 



d-l 

Note that the even part has a very simple form 
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Let us define h%[v ■ k,v 
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9aph 3 (v -k,v-q) + v a vp(- ■ •) 



or 



hs(v ■ k,v ■ q) = [h(v ■ k — v ■ q) — h(v 



(D. 



When the nucleon is on-shell, that is, v ■ k = v ■ q = 0, then h% becomes /13(C), 0) = — A(M 2 ). 
More generally, for small momentum, we have 

h 3 (v-k,vq) = -ACM 2 ) (2v-k—v-q) + -ri h(v ■ k) 2 - 3v ■ kv ■ q + (v ■ q) 2 }+- ■■ CD. 9) 

167T 3 L J 

where the ellipsis denotes finite and higher momentum terms. Finally consider h^Cv ■ q) 
defined by 

la 1/3 



ll (vl) 2 v(l-q) (l 2 -M 2 ) 
The /14 is a somewhat complicated function, 
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167r 3 

with h{y) given by ( p. 6 ) and the ellipsis again denotes finite and higher momentum terms. 
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Appendix E: Integrals for Two-Pion Exchange Currents 



Consider the integrals of the form 



Vs 



liv(l + k) (I 2 - M 2 ) [(I + q) 2 - M 2 ] 

which figure in two-pion exchange currents. For most of the cases, we do not need the terms 
proportional to as they appear multiplied by the spin operator and vanish. To utilize 
this, we assume that the spin operator is multiplied to the numerator. Now we have 
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where we have neglected terms proportional to v a or v 13 . After some algebra, we can get 
the following relations, 

q 2 Bi(k, q) + v ■ qB 2 (k, q) = h(v ■ k) — h(v ■ k — v ■ q), 
B 2 (k,q) = fi(q 2 )+v- q [B {k,q) - - 2v ■ kB (k,q). 

When v ■ k = v ■ q = 0, they become elementary functions, 

B (q 2 ) = -j^j o 1 dz^M 2 -z(l-z)q 2 , 
Bi(q 2 ) = 0, 

B 2 (q 2 ) = h{q 2 ) (E.3) 
For small but nonzero momentum, they become 
B (k,q) = -— + (v ■ k - ■ q^j rj H , 
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For two pion exchange graphs, we need to evaluate 
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(E.4) 
(E.5) 



v ■ I v ■ (I + q) (I 2 - M 2 ) [(/ + q) 2 - M 2 } ' 

In evaluating this function, we neglect terms proportional to v u , v a or vp because they 
vanish when multiplied by S u S a S 13 . With the parametrization explained in the text, we 
have, in the limit of v ■ q = 0, 
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where K\(q ) is defined at ( [B.5D , 



o 



For completeness, we also list results for vector currents for which we need the fol- 
lowing integrals 

(l + q )«lF>(2l + qy 



llvlv(l + q)(P- M 2 ) [{I + q) 2 - M 2 ] ' 
We first look at its low momentum behavior, 
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In the limit of v ■ q = 0, we have 
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Again we dropped terms proportional to v a or . 

Appendix F: Three-point Vertices (Figure 4) 

In this section as well as in the next, we classify graphs into Class A, Class V and 
Class AV. The graphs in Class V appear only with the vector current while the graphs 
in Class AV appear both for the vector and axial-vector currents. The graphs in Class 
A involving the axial-vector current do not figure in three-point vertices. We define two 
operators for the graphs in Class AV: = and T2 = 2gAS^. We write the expressions 
only for the axial-vector current for graphs in Class AV. The expressions for the vector 
current is obtained by interchanging and T% . Fig. 4c and Ad vanish because they are 
proportional to v ■ S. 



Class AVF 2 



r^(«) = y^A(M 2 ), (F.l) 

*%n(P) = J^S a T£S a h 3 (vk,vq), (F.2) 

*Knn(a) = ^r a q.SA(M 2 ), (F.3) 

iK NN (b) = ±zl Ta q-S^h 3 (vk,vq), (F.4) 



#22 The figure label a,b,c... is given in the parenthesis. 
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Class V 
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[q^B l {k 1 q)+v^B 2 {k,q)] 



+ {q-S, S»} B 1 {k, q) + [q- S, S"] 5 (fc, g)} 



(F.5) 



(F.6) 



T Vnn(9) = i£abc{q c - qbY^—pr 
1 



1 VlTTT 



(h) 



F 2 



{qc-qbYfM)+q»(qi-ql)fM 



(F.7) 
(F.8) 



Here the index a labels the isospin of the photon with four-momentum q a , the indices b and 
c the isospin of the pions with their momenta q b and q c and the momentum conservation is 

q a + qb + q c = o. 



Appendix G: Four-point Vertices (Figure 5) 

Here we define c/ M = q% + q£. For other notations, see Appendix F. Figures 5(«) — (n) 
vanish because they are proportional to t> • S. Here we restrict ourselves to the case of 
on-shell nucleons, v ■ k = v ■ (q a + q b ) = 0. 
Class A 



■^fi,ab 
nA 



(a) 



Class AV 



1 



■CabcTcV^fliq 2 ), 



2F 3 

-iS ab {2q a + 6q b r 
2g\ 

^abcTc j-,^ S v 



3F 3 

S-S+2{q-S) 2 







dq 2 



S-S+2{q-Sf^ 3 

d 
dq 2 



Bo(q 2 ) 
fl(q 2 ) + [q-S,S»]B (q 2 ] 



(G.l) 



(G.2) 



-^fj,,ab 
nA 



(d) 



5 

~12 

A. 

2F 3 



e abc r c ^A{M 2 ), 



e abc T c S a T>;S a A(M z ), 



(G.3) 
(G.4) 



Ka\c + /) = 4^3 { 4 *<^ " ■ 9a ^0 0> ■ 9a) + e abcTc [a(M 2 ) - 2t; • q a h£(v ■ q a )} } , 



■^fj,,ab 
nA 



(g + h) 



(SiS ab + e abc T c )-^S a q b ■ S T 2 M S a h A (v ■ q a ) 

g 3 

+(-3id ab - e abc T c )-^S a T£ q b - SS a h 4 (-v ■ q a ), 



(G.5) 



(G.6) 
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Class V 



9A 
2F 3 ' 



v ■ q a < (-2i5ab ~ eabcTc) 



l-S{2l + q a r 



- (2i5 ab - e abc T c ) 



lvl(P- M 2 ) [(/ + q a ) 2 - M 2 ] 
l-S(2l-q a y 



2i5 n A 



ivl{l 2 - M 2 ) [{I - q a ) 2 - M 2 } 
(l + q a )-Sg b -Sl-S(2l + q a y 
F 3 Ji v ■ (I + q a )v ■ I {I 2 - M 2 ) [{I + q a ) 2 - M 2 ] 



where 



B (q 2 ) = -^J Q 1 dz^M 2 -z(l-z)q 2 , 
h$(vq) = X - [ho(v ■ q) + ho(-v ■ q)} = ~yj 1 M 2 - (v ■ q) 2 , 
ho(v-q) = - [h (v ■ q) - h (-v ■ q)] . 



(G.7) 
(G.8) 

(G.9) 
(G.10) 
(G.ll) 



Now we study low-momentum expansion for on-shell nucleons with v • k = v • (q a + q b ) = 0. 
To second order in external momentum; 
Class A 



1 



2F 3 



EabcTc V 



a(m 2 ) - \e 



-iSab(q a +3q b r^M' 



(-abcTc 
+ ••• . 



2 A{ 

F 3 I 



v^S-SA(M 2 



S ■ Sq 2 + 2(q- S) 2 - [q-S,Sl*]M' 



Class AV 



■^fi,ab I 



kA ( c ) 

\d) 

Kt(e + f) 
Ka\9 + h) 



Class V 



-^e abc T c ^A(M 2 ), 
^e abc T c S a T?S a A(M 2 ), 



rpfl 

{-8i8 ab v-q a M' + e ah cTc [A(M 2 ) - 4(v • q a ) 2 rj\ } ■ ■ ■ , 



4i?3 

(-3iS ab + e ab cTc)^S a q b ■ ST£ S a 
+ {-Zi5ab ~ e ab cTc)^S a T$ q b -SS c 



M' + -vq a rj 

' . 2 

M - -v -q a V 



+ 



F»f(o +p + q) = ^v-q a [4i5 ab M'S» + e abc T c (S»v ■ q a + v*q a ■ S)rj\ +■■■, 

Kv( r ) = -2i5 ab ^{-2v^M'S a q b -SS a -(q a -Sq b -SS^-S^q b -Sq a -S) V } 
+ ■■■. 



We have used the notations M' = and rj = (i^jr) F(e). 



47 



Appendix H: Fermi-Gas Model for Two-Body Axial-Charge 
Operator 

Let \F) be the ground state of Fermi-gas model whose fermi-momentum is p? and 
\ph) = b p bh\F) be the one-particle (labeled by p) one- hole (labeled by h) excited state, 
where 6 a (6^) is the annihilation(creation) operator of a fermion state characterized by a. 
Consider the matrix element {ph\M.\F) (or its effective one body operator M e ff), 

(ph\M\F) = ( P \M cS \h) = J2 —<p,P\M\h,l3) (H.l) 

where gp is defined by |b a ,6^| = gp5 a fi. In computing this, it is convenient to define the 
antisymmetrized wave function \a,(3) in terms of the simple two-particle state \a,(3) 

so the matrix element (p, (3\M\h, (3) is of the form 

(p,P\M\h,0) = (p,0\M\h,P) - (p,(3\M\P,h). (H.3) 

The first term is the Hartree term and the second the Fock term. Rewriting \a) as \panT- a t a ) 
where p a is the momentum of the state labeled by a and m a (t a ) the third component of 
the spin (isospin) of the state a, we may write the axial charge operator as 

(l',2'|M|l,2) = ftmi.tXal [t^Mq) + T {2) Mq)} \hm u t 2 m 2 ) (H.4) 

where 7~W = i t\ X t 2 (cti + <r 2 ) • q, =i{f 1 + r 2 )(7i x cr 2 • q and q = p' 2 - p 2 = pi - p' 1; 

q = |q|- 

It is trivial to see that the Hartree term must vanish. Thus we are to calculate the 
Fock term. First we shall show that the matrix element of is equal to that of 
when summed over spin and isospin of occupied states. Doing the sum, we get 

H("Vp> m /3*/3|^ (1) \™>ptp, m h th) = -4 {m p t p \f (g> o\m h t h ) ■ q. (H.5) 



(p\M eS \h) = 4r ph a ph f (p - p„) 0(|p - p„|) (H.6) 



For T^ 2 \ we simply interchange the spin and isospin operators and get the same result. It 
then follows that the effective one-body operator of the axial-charge operator becomes 

dp 

/|p|<£f (2tt) s 

where <p{q) = 4>\(q) + <^ 2 (<?), Pe stands for the fermi-momentum of the ground state \F) and 
T p h(o~ph) is (t p \f\th) ((m p \ a \rrih)) ■ This form is particularly useful when the particle is on 
the fermi surface | p^ | = Pf, 

I 7%(P- Ph) H\P ~ PhD = -T^Srf f 1 dx (x 3 - x 5 ) <K2x ft ). (H.7) 

Ap\<Pf ( 2yr ) 47r JO 
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In order to give a meaning to this expression, we have to account for short-range correlations. 
Otherwise we can get erroneous results. For instance, a constant in momentum space (say, 
<j)(q) = constant) gives a contribution whereas it should be suppressed in reality. One way 
to assure a correct behavior at short-distance is to subtract the constant as one does for 
the Lorentz-Lorenz effect in 7r-nuclear scattering. However this procedure is not always 
practicable if one is dealing with non-polynomial terms. It is therefore preferable to go to 
coordinate space by Fourier-transform. For this, define f(r) by 

/(r) = /(^ eiq ' r0(9) - (K8) 

Using 

J lpl<PF (2vr)3 " 2^r Jl(PFrh 3l[X) ~ 

and J dQfe tp ' r = i4irpji(pr), 

we obtain 

(ph\M\F) = r ph a ph ■ p h / drrn^hiphr) —f(r). (H.9) 

7r Jo ar 

Introducing a correlation function g(r, d) where d is a parameter of g, we get the final 
expression 

(ph\M\F) = r ph a ph ■i>h J ^ L J d drrj l { Pe r)j 1 { Ph r) -f(r). (H.10) 

In the numerical results discussed in the next, we have used the simplest correlation function, 
g(r,d)=9(r-d). 



Appendix I: The Role of Vector Mesons 

In this Appendix, we describe briefly the role that vector mesons play in the axial- 
charge transitions in chiral perturbation theory. We will in particular establish that vector 
mesons can contribute only at O(^) and hence their contributions are suppressed to the 
chiral order we are concerned with. For simplicity, we shall consider the vector field only. 
The axial vector field a±n could also be included but we shall leave it out since it plays even 
less significant role in our case. Let 

V„ = t a V°, Ti(t a t b ) = h ab (1.1) 

denote the spin-1 field. The index a and b are (1, 2, 3) for SU(2) with t = J, and (0, 1, 2, 3) 
for U(2) with to = \- The a = (1, 2, 3) components correspond to the p mesons and a = 
to the u) meson. We write the relevant part of the Lagrangian as* 23 

C = N[ 1 »(id„ + gV„ + g Al5 iA„)-m}N+^(iA^iA») + ±M 2 F 2 (Z) 

* 23 We have not included the axial- vector field ai, although it is not difficult to do so. For the axial-charge 
process we are considering the ai field does not play an important role. For the Gamow- Teller operator, 
however, the axial field may not be ignorable. 
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+ \m$((v» - ~r M ) ) - \{v, u vn + c an (i.2) 

where 

V^ = d^V u -d v V^-ig[V^V u ], (1.3) 
and r„ (A„) were given in Section 3 and explicitly take the form 

iT^ = t- 
iA^ = t- 

where (An) is the external vector (axial- vector) field and the ellipsis denotes terms 
involving more than three fields. Here C an is an "anomalous parity" piece involving the 
totally antisymmetric e tensor which we do not explicit here as it does not contribute. Also 
four-fermion interaction terms do not figure in the discussion. In (|L2}), the constants g 
and My can be identified as the VNN coupling constant and the mass of the V meson 
respectively. We are using the short-hand notation 

(X)^2Tr(X) (1.5) 

for any X. This convention is convenient due to the normalization of t a , (XY) = X a Y a for 
any X = t a X a and Y = t a Y a . 

Before proceeding, let us note a few characteristics of this Lagrangian: 

• It is vector gauge-invariant (or hidden gauge invariant) provided that gV^ transforms 
as iT^ does, — > UV^U^ — ^d^U • UK It is also invariant under (global) chiral 
transformation apart from the pion mass term. 

• It has vector-meson dominance. 

• When My goes to infinity, we recover our previous chiral Lagrangian involving only 
7r's and nucleons N. 

• There is a V"f mixing but the mixing is trivial in the sense that the photon field 
appears only as an external (non-propagating) field. 

For the reasons spelled out in the main text, we wish to transform the Lagrangian to 
a form appropriate for heavy fermion formalism. Including the "1/m" terms, we have 

F 2 1 
C = B(ivD + 2g A S-iA)B + —(iA^iA^) + -M 2 F 2 (^) 

+ J-B (-D 2 + (v ■ D) 2 + [5", S v ] D v ] - g\(v ■ iA) 2 - 2ig A {v ■ iA, S ■ D}) B 
+ l -MU(v,- l -V^)-\{V llu V^) + C, an (1.6) 



V M + — t? x An 
b 



^7? X (9 U 7T + 



1 - 1 

An + — 7T X Vfj, - + 



(1.4) 
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where = — igV^. Now let us calculate the tree-order contribution of the vector 
mesons to the two-body axial charge operator. Three types of graphs contribute. The 
relevant graphs are given in Figure 11. 

First we find that the graph (c) does not contribute. To see this, note that G-parity 
does not allow the couplings ppA^ and ujloA^. The coupling for A^puj is of the form 



tPp coming from the anomalous-parity term C an of (|L6|). In the figure (c), each 



vector meson brings in as one can see in (1.6), so that we have e^ vol ^v a vp = 0. 
Working out the graphs (a) and (6) f^ 24 \ we get 



^(1) 



A^{2) 



%T\ X T 2 



9A 



Ml 



1 



2F 2 M 2 



2 M 2 



1 



(jl + T 2j 



m N 
9A 



Qi 

[S?,qi ■ Si 



4 



<?2 ■ S 2 



v ■ Q2 

m N 



vi ■ gi n 



Ml 



2F 2 Ml 



q 2 ^ m N 



+ (1^2) 



+ (l<->2) 



(Ii 



v 2 ■ S 2 H [Si ■ S 2 ,q 2 ■ S 2 ] 



m N 



v 2 -q 2 



My 



v 



(1.9) 



where we have used the KSRF relation My 



2g 2 F 2 and defined 



m N 



(1.10) 



and qi 



P'i 



have (setting q 2 - 

A°(l) = %T\ x T 2 



Pi, i = 1,2. Now noting that v ■ q 
-<7i = q) 



9A 



Ml 



2F 2 Ml 



+ (1 



A°(2) 



(n + T 2/ 



9A 
2F 2 



O 



Q 2 



m 



N , 



1 



+ o 



and S? = O 



M 2 \M 2 - q 2 



Ml 



Q 2 



m 2 N Ml 



q ■ S 2 



l + O 



mjv 



Q 2 



m 



N , 



we 



(1.11) 
(1.12) 



The leading part of Equation fll.ll ) is nothing but the one-pion exchange current with one- 
loop radiative corrections ( |9l| ) expressed now in terms of a vector-dominated Dirac form 
factor FY ■ (In fact, with the Lagrangian (|I.6|), the soft-pion contribution corresponds to 
fll-liP in the limit My — > oo.) There is no further correction to what has already been 
obtained with our Lagrangian given in its full glory in Appendix A. This corroborates our 



#24 



In the spin-1 propagator 



q 2 - AL 



9 + Ml 



(1.7) 



the term proportional to 773- <? q is a correction term of O ( 773- ] relative to the leading term (oc p M ). It 

m v \ M v I 

is further suppressed in the vector-meson exchange between two nucleons since the VNN vertex function is 

proportional to u M and v ■ q = O (^%T^j ■ ^° effectively this term is of order J^ M 'i and hence can be dropped. 
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(1 2) 

argument that the counter terms k\ ' cannot come from one vector-meson exchange in the 
limit my — * oo. This also establishes our assertion that vector mesons do not modify our 
result on the "chiral filter mechanism." 

For completeness, we give the corresponding axial-charge operator in coordinate space 



Ml^ + M\r = (1 + Ssoft) ± [(m + i) e~ Mr - (m v + I) 



e -Myr 



(1.13) 



where 5 so f t = M -I^_ M 2 an d we have dropped the terms of order m N 2 . This is the vector- 



dominated form of Aii-n-, in place of ( |104| ): The second term of vector-meson range in ( |1.13 ) 
is the counterpart to the shorter-ranged loop correction in ( |104D . In fermi-gas model fll.13 ) 
predicts roughly the same quenching as the loop calculation (|104|). 



52 



FIGURE CAPTIONS 



Figure 1 

Generic nuclear electroweak currents up to two body. The solid line represents 
the nucleon, the blob with a cross the coupling of electroweak fields and the 
shaded blob without cross stands for the strong interactions. 



Figure 2 

Two-body exchange currents: (a) One-pion exchange; (b) two-pion exchange. 
The solid blob represents a strong-interaction vertex and the shaded blob with 
a cross the vertex involving an external field and strong interactions. The solid 
line represents the nucleon and broken line the pion. 



Figure 3 

One- loop graphs contributing to the nucleon self-energy S. As in Figure 2, the 
solid line represents the nucleon, the broken line the pion. 



Figure 4 



One-loop graphs contributing to the three-point G^NN vertex where G^ = 
(Vu) is the external axial-vector (vector) field, the encircled cross representing 
the field coupling. Here and in Fig. 5, vector- field couplings are also drawn for 
comparison and for later use in [14]. 



Figure 5 

One-loop graphs contributing to the four-point G^ttNN vertex. For axial-charge 
transitions, only the graphs (a)-(f) contribute. 



Figure 6 

One- loop graphs contributing to two-body two-pion exchange currents ((a) — 
(h)), four-fermion-field contact interaction currents ((i) — (j)) and "recoil" cur- 
rent (k). The pion propagator appearing in (a) — (jj) is the Feynman one while 
that in (k) is a time-ordered one. Only the graphs (a), (b), (c) and (d) survive 
for the axial-charge operator. 
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Figure 7 



47rr 2 A4f ree (solid line) and Anr 2 M.i oop (broken line) denned in Eqs.(|l0 
(106) vs. r in fm. Here and in Fig. 8, we have set = = i. 



and 



Figure 8 

r [ji(PFf)] 2 M-tree (solid line) and r[ji(pFr)] 2 M.[ oop (broken line) vs. r withpi? 
1.36 fm^ 1 (corresponding to nuclear matter density). See the caption for Fig. 
7. 



Figure 9 

The ratios of the matrix elements n^ t x ^ i n fermi-gas model vs. p/po for d = 
0.5, 0.7 fm for X = lir, 2tt, lir + 2ir corresponding to one-loop correction to the 
one-pion exchange graph, one-loop two-pion exchange graph and the sum of the 
two, respectively. 



Figure 10 

Three-body currents: a) Genuine three-body current with Feynman pion propa- 
gators; b) "recoil" three-body currents with time-ordered pion propagators; the 
ellipsis stands for other time-orderings and permutations. Both (a) and (b) are 
of order 0(Q 3 ) relative to the leading soft-pion term. 



Figure 11 

Vector-meson contribution with the Lagrangian (|Lq) to the two-body axial 

charge operator. V and V stand for vector mesons of mass My. For the 
axial current, V = p and V = oj. 
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